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ABSTRACT .  A  literature  survey  and  an  evaluation  ot  applicable 
literature  concerning  pressure  and  heat-transfer  distributions  over 
bodies  of  revolution  in  high-subsonic  to  low-supersonic  flight  is 
presented.  Methods  for  determining  the  pressure  distribution  on 
bodies  of  revolution  are  discussed  and  compared  with  experimental 
data.  Variations  of  the  modified  Newtonian  method  may  be  applied 
to  a  large  class  of  bodies  of  revolution  for  free-stream  Mach 
numbers  greater  than  0.75. 

Appropriate  experimental  laminar  and  turbulent  heat-transfer 
data  for  a  blunt-body  vehicle  at  flight  Mach  numbers  between  1.6  and 
5  are  compared  with  calculated  values.  Experimental  data  could  not 
be  found  for  blunt-body  turbulent  heat  transfer  in  the  free-stream 
Mach  number  range  from  0.7  to  1.6. 

An  extensive  bibliography  covering  prediction  techniques  and 
experimental  data  on  pressure  and  heat-transfer  distributions  is 
included  as  the  appendix. 
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NOTATION 

Velocity  gradient  at  the  stagnation  point 

Local  skin-friction  coefficient,  2t  /p  u  2 

wee 

Local  Stanton  number,  q/p  u  [h  -  h  ) 

T*  e  e  \  aw  w) 

Local  Stanton  number  based  on  reference  conditions, 

■Vc\{haw  -  h„) 

Pressure  coefficient  =  ^P  -  Pocj/(y/2)  P^M^2 

Specific  heat  at  constant  pressure 

Body  diameter  at  base  of  nose  section 
A  function,  defined  by  context 
Enthalpy 
h/h 

s 

e 

Adiabatic  wall  enthalpy 

Mjd/L) 

Thermal  conductivity 
Length  of  nose  section 
Length  of  body 
Mach  number 

Local  Nusselt  number  based  on  reference  conditions, 

<L  U  )  s/k  /h  -  h  \ 
w  \  p  7  \aw  w/ 

Pressure  at  the  surface  of  the  body 

Prandtl  number,  0.71,  yc  /k 
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Heat  transfer  per  unit  area  per  unit  time 


Radius  of  body  cross  section  normal  to  the  axis  of  symmetry 


Re  0 


h 


Reynolds  number,  p  u  R  /y 
s  00  o  s 
e  e 


Local  Reynolds  number  based  on  reference  conditions, 
p  u  s/y 


Hemisphere  radius 


r  /R 
o  o 


Body  radius  at  base  of  nose  section 


Radius  of  blunting  at  tip  of  nose  section 
Radial  distance  from  axis  of  symmetry 


Area  of  cross  section  normal  to  axis  of  symmetry 

Wetted  length  from  stagnation  point 

s 

R 

o 

Temperature 

Velocity 

u/u 

e 

Perturbation  velocity,  x  component 

Perturbation  velocity,  y  component 

Perturbation  velocity ,  z  component 
Axial  distance  from  nose 

i 

I 

1 
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y  Ratio  of  specific  heats 

6  Angle  between  flight  direction  and  radius  vector  from 

center  curvature  of  nose 

0^  Local  slope  of  body  surface  relative  to  the  axis  of 
symmetry  (the  angle  between  tangent  to  body  surface 
and  the  axis  of  symmetry) 

y  Viscosity 


v 

P 

P 


T 

T 


Prandtl-Meyer  angle 
Mass  density 


P  /Pc 
e  s 


w 


P  /p  h 
s  s  s 
e  e 


Characteristic  body  parameter,  d/£ 
Wall  shearing  stress 


$  Velocity  potential 

<t>  First-order  perturbation  velocity  potential 

0  Second-order  perturbation  velocity  potential 
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Subscripts 

e  Conditions  external  to  the  boundary  layer 
N  Conditions  at  nose  of  pointed  body 

a  Stagnation  conditions 

t  Total  conditions 

v  Wall  conditions 

00  Free-stream  conditions 

2  Conditions  evaluated  behind  a  shock  wave 

Superscripts 

* 

Denotes  conditions  evaluated  at  the  state  corresponding  to  the 
reference  enthalpy 
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Section  1. 
INTRODUCTION 


The  purpose  of  this  report  is  to  present  a  review  of  the  available 
methods  for  determining  the  pressure  distribution  on  bodies  of  revolution 
in  the  subsonic,  transonic,  and  supersonic  flight  regimes  and  a  comparison 
of  the  results  of  these  methods  with  available  experimental  pressure 
distributions.  In  Section  2,  the  methods  discussed  are  those  that  are 
most  useful  as  "engineering"  tools  and  provide  the  pressure  distribution 
with  a  minimum  of  effort.  These  methods  are 

1.  Modified  Newtonian  theory: 

a.  An  empirically  modified  approach  applicable  in  the 
range  M  *  0.75  to  M  *2.0 

b.  A  hybrid  method  combining  the  modified  Newtonian  and 
shock-expansion  theories  for  the  range  M^  >  1.5 

2.  Perturbation  theory  as  developed  by  Van  Dyke  and  Sprieter 

3.  Similarity  rules. 

Methods  2  and  3  may  be  applied  over  the  entire  range  of  Mach  numbers. 
However,  these  methods  require  considerably  more  effort  than  does 
Method  1. 

In  addition,  brief  reference  is  made  to  the  most  successful  of 
computer  techniques  used  for  determining  pressure  distribution  for 
blunt  bodies  in  supersonic  flow. 

Section  3  presents  a  review  of  the  engineering  methods  of  analysis 
for  the  convective  heating  of  blunt-body  vehicles  in  the  high-subsonic 
to  low-supersonic  flight  regimes.  Comparison  of  these  simplified 
analytical  methods  with  appropriate  experimental  data  its  also  presented. 

The  objective  of  Section  3  is  to  present  methods  of  prediction  that 
are  both  convenient  in  their  utilization  and  accurate  in  their  agreement 
with  available  experimental  data.  These  methods  of  prediction  are  (1) 
of  convective  heat  transfer  at  the  stagnation  point  of  a  blunt-body 
vehicle,  and  (2)  of  distribution  of  the  heat  transfer  over  the  surface 
of  the  vehicle.  With  this  objective  in  mind,  a  literature  search  was 
performed  to  locate  both  the  methods  of  analysis  and  the  appropriate 
experimental  data. 


1 


NWC  TP  4570 


Considerable  effort  has  been  devoted  to  determining  the  pressure 
distribution  over  bodies  of  revolution  in  subsonic  and  supersonic  flight. 
With  the  exception  of  a  few  notable  early  efforts,  the  bulk  of  the  work 
has  been  done  during  the  past  20  years.  The  appendix  of  this  report 
presents  a  bibliography  covering  both  theoretical  developments  and 
available  experimental  pressure  distribution  data.  This  bibliography 
comprises  567  entries. 
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Section  2 


PRESSURE  DISTRIBUTION  ON  BODIES  OF  REVOLUTION 


MODIFIED  NEWTONIAN  THEORY 

The  central  hypothesis  of  Newtonian  flow  theory  is  that  the  component 
of  momentum  of  the  oncoming  flow  normal  to  the  body  surface  is  transferred 
to  the  body  while  the  tangential  component  remains  unchanged.  Various 
authors  (Ref.  1)  have  pointed  out  that  Newtonian  theory  is  reached  only 
through  a  limiting  process  as  Mo  +  “  and  y  1.  Cole  (Ref.  2)  presents 
a  more  detailed  examination  of  the  relationship  between  Newtonian  theory 
and  gas  dynamics;  but  again,  he  concludes  that  Newtonian  theory  is 
applicable  only  for  large  Mach  numbers.  Based  on  these  conclusions, 
Newtonian  theory  was  considered  only  as  a  limiting  case  and  was  originally 
used  only  to  describe  pressure  distributions  in  the  hypersonic  flight 
regime  (Ref.  1,  2,  3).  However,  it  was  soon  found  from  experimental 
data  that  in  a  modified  form,  application  was  valid  also  In  the  supersonic 
flight  regime  (Ref.  4-9). 


In  its  original  form  the  Newtonian  concept  provided  a  coefficient 
of  pressure  in  the  form 


2  sin20. 


This  relationship  indicates  dependence  only  on  the  local  slope  of  the 
surface,  0b-  Agreement  with  experiment  was  not  achieved,  however,  until 
Lees  (Ref.  3)  suggested  a  very  useful  modification  of  this  expression  in 
the  form 


sin20. 


In  Eq.  2,  C  is  the  maximum  value;  of  the  coefficient  cf  pressure 
^max 

on  the  body  surface.  For  a  blunt  body  in  supersonic  flow,  the  maximum 
value  of  Cp  is  at  the  stagnation  point,  and  its  value  may  be  conveniently 
computed  from  the  following  expression: 


.  (P'2/P")-1 

<J  =  ~ 

*max  (y/2)  M  2 
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where  P  / P  is  determined  from  the  Rayleigh  formula: 
2 


P.  /Poo 
2 


,  Y/Y-l 


’(y+DmJ  ' 

1  •  1  — 

1 

Y  +  1 

2 

2YMra2  -  (Y-l)  . 

1/Y-l 


(4) 


Values  of  C  computed  by  using  Eq. 


3  and  4  are  compared  with  experi- 

1  max 

mental  values  in  Fig.  1.  Figure  1  demonstrates  an  excellent  agreement 
between  theory  and  experiment.  As  becomes  large,  computed  values  of 
C  for  air  approach  1.83  rather  than  2. 

^max 


FIG.  1.  Maximum  Coefficient  of  Pressure  for  Blunt  Bodies 
of  Revolution  as  a  Function  of  Mach  Number. 


Oliver  (Ref.  10)  presents  the  following  simple  formula  which  may  be 
used  for  determining  C  if  Ma  >  2: 

^max 
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C 

p 

max 


X±1  L  _  _2_  JL. 

Y+1  Y+3  Mro2 


(5) 


From  Eq.  2,  3,  4,  and  5,  Cp  is  seen  to  be  a  function  of  Mro  and  y,  as  well 
as  of  %.  Comparison  of  Eq.  2  with  experimental  data  for  hemispherically 
capped  cylinders  is  presented  in  Fig.  2.  The  Mach  numbers  of  the  data 
ranged  from  1.97  to  4.76.  Equation  2  represents  the  data  quite  adequately. 


LOCAL  SLOPE,  6b  OEG 


FIG.  2.  Distribution  of  Pressure  Coefficient 
Ratio  Ovei  a  Hemispherical  Nose,  Mro  >_  1.97. 


For  bodies  of  revolution  which  are  only  slightly  blunted  or  for  a 
pointed  body,  the  modified  Newtonian  method  may  also  be  applied  if  C 
is  computed  as  suggested  by  Love  (Ref.  6):  Pmax 


max 


C  / sin2  0, 

pN/ 


(6) 
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where 


C 

P 


from  tables  of  conical  flow  (Ref.  11)  for  a  cone  semivertex  angle,  0jq, 
corresponding  to  the  local  slope  at  the  nose  of  the  pointed  body. 


An  example  of  what  is  meant  by  an  "only  slightly  blunted"  body 
would  be  a  hemispherically  capped  ogive  with  a  cap  radius  that  is  small 
relative  to  the  maximum  cross-section  radius  of  the  ogive  (say  rc/r^  less 
than  approximately  0.2).  For  such  a  case,  the  Cp  at  the  stagnation  point 
will  correspond  to  that  of  Eq.  3,  but  the  distribution  away  from  the 
stagnation  point  must  be  determined  on  the  basis  of  the  C  given  by 

^max 

Eq.  6  for  an  equivalent  pointed  body  (Ref.  12).  By  way  of  explanation, 
an  equivalent  pointed  body  is  that  pointed  body  which  would  exist  if  the 
nose  were  not  blunted. 


Love  (Ref.  6)  presents  a  very  interesting  comparison  of  this 
generalized  Newtonian  theory  (Eq.  2  and  6)  with  experimental  data  for 
a  variety  of  bodies  of  revolution  over  Mach  number  ranges  from  2  to  7. 
This  comparison  is  reproduced  in  Fig.  3  with  the  scatter  in  the  data 
Indicated  by  the  band  width. 

The  modified  or  generalized  Newtonian  prediction  begins  to  deviate 
from  experimental  values  as  the  local  slope  approaches  zero,  or  as  a 
point  of  discontinuity  in  slope  is  approached.  Lees  and  Kubota  (Ref.  13) 
suggest  that  in  such  a  region,  a  local  shock-expansion  method  be  used  to 
predict  the  distribution.  Tc  accomplish  this  the  modified  Newtonian 
equation 


C 

P 


C  pin2  0v 

P  D 

max 


6 


1.0  0.9  0.8  0.7  0.6  0.5  0.4  0.3  0.2  0.!  0.0 

SIN2  8. /SIN2  8 
d  n 


FIG.  3.  Love's  Correlation  for  Generalized 
Newtonian  Theory. 

j 


t 


and  the  Prandtl-Meyer  relation 

p  (d!0  *  -  1 


(7) 


are  matched  at  the  point  on  the  nose  surface  where  both  the  pressure 
and  pressure  gradient,  predicted  by  these  two  solutions,  are  equal. 
This  method  was  further  investigated  by  Wagner  (Ref.  14),  Jones 
(Ref .  15) ,  and  Bllck  and  Francis  (Ref .  16) .  Blick  and  Francis  suggest 
the  following  formula  for  determining  the  matching  angle  in  the  Mach 
number  range  2  to  15: 


9 

m 


0.349  + 


198  +  0.80 


(Mro  -  2.8) 


-0.455 


-1 

cos 


1.05/ (In 


(8) 
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Vendemla  (P.ef.  9)  presents  an  excellent  summary  of  this  hybrid 
modified  Nevtonian/shock-expansion  method.  A  number  of  comparisons 
of  the  method  with  experimental  data  or  "exact"  theoretical  methods 
for  ogives,  Von  Karman  minimum  drag,  and  3/A  power  law  bodies  in  the 
Mach  range  *  1.5  to  6  and  for  hemispheres  in  the  =  1.82  to  A. 76 

range  are  presented  in  Ref.  9.  Two  of  these  comparisons  are  reproduced 
as  Fig.  A  and  5.  These  figures  show  that  the  agreement  between  the 
present  theory  and  experiment  or  the  more  "exact"  theory  is  excellent. 


FIG.  A.  Comparison  of  Modified  Newtonian/Shock- 
Expanslon  Method  With  Method  of  Characteristic 
Solution  for  Tangent  Ogive  Bodies  at  Supersonic 
Mach  Numbers . 


Vendemla  also  developed  the  following  expression  to  determine  the 
matching  point  for  the  modified  Newtonian  and  shock-expansion  methods. 
From 


C 

P 


C  sin2  9, 
p  b 

max 


<Y/2>PcoMoo2 


8 
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BODY  PROFILE:  r  -  Vi  -  7/2  SIN2  + 

V* 


0  0.1  0.2  0.3  O.R  0.5  0.6  0.7  0.8  0.9  1.0 


NONDIHENSIONAL  AXIAL  DISTANCE,  x/L 

FIG.  5.  Comparison  of  Modified  Newtonian/ 
Shock- Expans  ion  Method  With  Experimental 
Data  for  Von  Karman  Minimum  Drag  Body. 


it  follows  that 


1  _dP 

p  de, 

D 


(y/2)  M  2C  sin  26. 

00  p  b 

_ _ _ max _ 

1  +  (y/2)  M  2  C  sin2 6, 

00  p  b 

max 


(9) 


Computations  are  then  made  to  determine  the  matching  point  0m  where  P 
and  dP/d0b  are  the  same  from  both  equations.  Vendemia  then  presented 
the  results  in  graphical  form  as  a  plot  of  the  matching  angle  0m  versus 
M^  for  hemispheres  and  hemisphere  cones.  The  result  is  reproduced  in 
Fig.  6.  It  should  be  noted  that  for  low  free-stream  Mach  numbers,  there 
is  no  point  at  which  both  the  pressure  and  pressure  gradient  match,  for 
the  two  methods  and  the  curve  of  Fig.  6  have  been  extrapolated  in  this 
area.  Vendemia  (Ref.  9)  found  that  the  matching  point  on  a  tangent 
ogive  can  be  determined  from  the  following  empirical  equation: 

(x7L')m  =  0.60  +  0.15K  for  K  <  2.667  (10) 


where 

K  -  (d/L) 

9 
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This  equation  may  also  be  applied  to  slightly  blunted  ogives  by  relating 
the  (x'/L'^  computed  from  Eq.  10,  which  would  be  the  nondlmenslonal 
distance  from  the  noue  of  the  equivalent  pointed  ogive  to  the  correspond 
lng  point  on  th*  blunted  ogive. 


For  blunted  ogives  where  rc/r^  is  considerably  greater  than  0.2, 
an  alternative  method  of  determining  the  pressure  distribution,  to  be 
discussed  shortly,  will  be  more  effective. 


FIG.  6.  Matching-Point  Angle  for  Hemisphere 
or  Hemisphere-Cone  Body.  (Vendemia  graph. 
Ref.  9.) 


For  predicting  the  pressure  over  a  cylindrical  afterbody,  Vendemia 
(Ref.  9)  gives  the  following  equation: 


10 
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C 

P 


C 

P(x/L-1) 


(ID 


where  Ax/L  is  the  nondimensional  distance  from  the  base  of  the  nose 
where 


x/L  =  1 

The  method  outlined  above,  using  Eq.  2,  6,  10,  and  11,  has  been 
applied  to  a  hemispherlcally  capped  ogive  body  and  is  compared  with 
experimental  data  (Ref.  12)  in  Fig.  7  and  8.  These  figures  illustrate 
the  usefulness  of  this  method  in  obtaining  good  results.  It  should 
also  be  noted  that  Eq.  3  again  gives  the  correct  value  of  Cp  at  the 
stagnation  pclnt.  fhe  calculations  for  Fig.  8  are  given  as  Example  1, 
page  17. 


BODY  PROFILE  ^  * 


FIG.  7.  Pressure  Distribution  on  a  Spherically  Tipped 
Tangent  Ogive  Eody. 


A  further  significant  application  of  modified  Newtonian  theory  was 
carried  out  by  Andrews  (Ref.  8).  For  the  case  of  a  hemisphere  cylinder, 
an  empirical  modifying  function  was  added  to  the  modified  Newtonian 
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FIG.  8.  Comparison  of  Modified  Newtonian/Shock-Expansion 
and  Van  Dyke's  Second-Order  Method  With  Experimental  Data 
for  a  Spherically  Tipped  Tangent  Ogive  Body. 


expression,  Eq.  2,  and  agreement  was  achieved  with  experimental  data  from 
Ref.  8  and  17  for  a  Mach  number  range  from  0.75  to  2.0.  The  equation 
used  was 


C 

P 


C  sin2 

p 

max 


eb  +  F(VK) 


(12) 


where 


^0b,M^  -  |o.78  M^-2* 27  sin  eb  -  0.95  exp[-2.235  (1^-1)]}  cos  ^ 


This  equation  Is  presented  in  Fig.  9  and  10  and  compared  with  experimental 
results  (Ref.  'S,  8,  9,  and  17). 

In  the  subsonic  regime,  C  for  blunt  bodies  is  computed  from  the 
isentropic  relation  Pmax 


c 

Pmax  (y/2)  M^2 

This  equation  is  compared  with  experimental  results  in  Fig.  1. 


(13) 
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LOCAL  SURFACE  SLOPE,  6.  DEG 

D 

FIG.  9.  Comparison  of  Eq.  12  With  Experimental 
Data  for  Hemisphere  Cylinders. 


A  modified  form  of  Eq.  11  may  also  be  used  ix  combination  with 
Eq.  2  or  12  to  predict  the  pressure  distribution  on  a  hemisphere  cone 
combination  in  the  supersonic  regime.  The  pressure  distribution  for  the 
hemispherical  nose  is  determined  by  using  either  Eq.  2  or  12  with  Eq.  3. 
The  cone  tables  (Ref.  11)  are  then  used  to  determine  the  pressure  coeffi¬ 
cient  for  a  cone,  C  ,  with  a  semivertex  angle,  8^,  equal  to  the  slope 
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FIG.  10.  Comparison  of  Eq.  12  With  Experimental 
Data  for  Hemisphere  Cylinders. 


of  the  conical  section  of  the  body  of  interest.  This  will  be  the  value 
of  C  reached  on  the  cone-  surface  some  distance  downstream  of  the 
Pc 

hemisphere-cone  junction.  The  pressure  distribution  on  the  conical 
section  can  then  be  approximated: 


C  -  AC 
P  P(j-c) 


(14) 


14 


ft 


,  w 


JfWC  TP  4370 

o 


where 


AC  -  C  -  C 

P(j-c)  Pj  Pc 


=  Cp  from  Eq„  2  at  the  hemisphere-cone  junction 


C 

P 


c 


for  a  cone  of  0„ 


The  use  of  this  approach  is  illustrated  by  comparison  vrith  experimental 
data  (Ref.  18)  in  Fig.  11  and  12. 

For  a  blunt  body  with  other  than  a  hemispherical  nose,  Eq.  12  may 
be  used  with  good  results  if  the  function  F  is  multiplied  by  the  fineness 
parameter  (r/L): 

C  -  C  sin2  6fa  +  (rfe/L)  F(6b,  M)  (15) 

"  Pmax 


r^  *  radius  at  the  base  of  the  nose  section 


L  =  axial  distance  from  stagnation  point 
to  nose  base 


F(8, ,  M)  =  same  as  that  in  Eq.  12 
b 


Eq.  15  indicates  that  tv  e  influence  of  the  modifying  function  F(t>b,  M^) 
decreases  as  the  body  becomes  more  slender. 

This  equation  was  used  tL  compute  the  pressure  distribution  on  a 
body  having  an  elliptical  nose  with  r^/l.  -  U.61.  The  results  are 
compared  with  experimental  data  (Ref.  19)  in  Fig.  13.  Equation  15  was 
also  used  to  predict  the  distribution  on  a  prolate  spheroid  (ellipsoid 
of  revolution)  with  r^/L  =  0.167,  and  the  results  are  compared  with 
experimental  data  (Ref.  20)  in  Fig.  14.  From  Fig.  14,  note  that  the 
pressure  gradient  is  predicted  quite  well,  even  though  the  actual 
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EXPT.  M^  REF. 

0  1-20  18 

D  ANDREWS  (EQ.  1  2  ) 

(EQ.  H) 

C  FROM  CONE  TABLES 
.  P  \ 

V 

o 

\  n  — JV-w 

\  0 

0.1  0.2  0.3  0.4  0.5  1.0  1-5  2.0  2.5  3.0 


DISTANCE  FROM  NOSE,  x  IN. 


FIG.  11.  Pressure  Distribution  on  a  Hemispherically 
Capped  15-Dey  Cone. 
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0  0.1  0.2  0.3  o.k  0.5  1.0  1.5  2.0  2.5 

i—  HEMISPHERE - 4- - -  CONE - 

DISTANCE  FROM  NOSE,  x  IN. 


FIG.  12.  Pressure  Distribution  on  a  Hemispherically 
Capped  15-Deg  Cone. 
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FIG.  13.  Pressure  Distribution  on  an  Ellipsoid  of 
Revolution-Nosed  Cylinder. 


magnitude  of  Cp  from  Eq.  15  and  experiments  differ  somewhat.  As  an 
example  of  this  approach,  the  calculations  corresponding  to  Fig.  13 
are  presented  on  page  27. 

The  chief  advantage  of  the  Newtonian  method  is  that  with  a  minimum 
of  effort,  good  accuracy  is  available. 

Table  1  presents  a  summary  of  the  available  pressure  distribution 

data. 


EXAMPLE  1:  MODIFIED  NEWTONIAN/SHOCK-EXPANSION  METHOD 

As  a  numerical  example  of  the  modified  Newtonian/shock-expansion 
method,  consider  again  the  data  shown  in  Fig.  7  for  M^  ■  2.17.  The  body 
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DISTANCE  FROM  NOSE,  t 


FIG.  14.  Comparison  of  Prediction  Techniques  With 
Experiment  Pressure  Distribution  on  a  Prolate 
Spheroid.  (Data  from  Ref.  24.) 


is  a  spherically  capped  ogive  with  rc/rj,  *  0.158.  By  assuming  that  this 
body  may  be  treated  in  terms  of  an  equivalent  pointed  body,  the  pressure 
distribution  may  be  computed  as  follows: 

1.  The  coefficient  of  the  pressure  at  the  stagnation  point  (Fig.  1) 
is  given  as  C  -  1.70. 

pc 

2.  To  determine  the  distribution  away  from  the  stagnation  region, 
first  determine  8^  for  the  equivalent  pointed  body.  For  an  ogive, 


sin  0 


L  -  x 


R 


where 


L  ■  length  of  ogive 
x  =  axial  distance  from  the  nose 

R£  ■  radius  of  ogive  surface  curvature  (in  the  x-y  plane) 
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TABLE  1.  Pressure-Distribution  Data 


Data 

Source 

Body 

Cooflouritlon 

Nidi 

Nurter 

- Tist - 

Reynolds 

Nuirter/pt 

Theoretical 

Comparison 

Andrews ,  J.  S. 

(AO  465-140,  Boeing  D2- 
22947-1;  1964) 

Hemisphere  cylinder 

0.8 

0.9 

1.0 

1.076 

5.5  •  10* 

6.0  •  10* 

6.1  10* 

6.2  •  10* 

Equation  12. 

(Imported  lo  And  rev*  but 
original  source  Is  Morris 
and  Grobll  (16)) 

1. 59 

1.79 

Not  reported 

Baer,  A.  L. 

(AO  261-501,  AEDC  TV  61- 
96;  1961) 

Heil sphere  cylinder 

1.99 

5.00 

4.03 

5.06 

6.03 

6.10 

2.02  -  10s 

4.2  •  10* 

5.3  •  10* 

6.1  •  10* 

4.4  •  10* 

2.5  •  10* 

Modified  Newtonian 
Equation 

Buglla,  J.  J. 

(NASA  TN  D955;  1961) 

Hemisphere  cone 
'coo.  •  l4'30’ 

Free  flight 
1.97 
to 

3.14 

Approximately 

10  ■  10* 

Modified  Newtonian 

Carter,  H.  S.  and  Breasette, 

V.  E. 

(NACA  1M  L57C18;  1957) 

Hemisphere  cone 

2.0 

14  •  10* 

None 

Flat- Faced  cone 

2.0 

14  •  10* 

Chauvln,  L.  T. 

(NACA  KM  L52K06;  1952) 

Hemisphere  cone 

2.05 

2.54 

3.04 

4.4  •  10* 

4.57  •  10* 
4.16  •  10* 

None 

Eatabrook.  B.  B. 

(AEDC  TV- 59-12,  AD  216 

698) 

20*  Cone  Cylinder 
and  parabolic  arc- 
cy  Under 

0.70 

0.80 

0.90 

0.95 

0.975 

1.00 

1.025 

1.050 

1.10 

1.15 

1.20 

1.30 

1.40 

Benge — from 
1.36  a  10* 
to 

4.0  x  10s 

Compared  with  the¬ 
oretical  cone  data 

Gardiner,  G.  K. 

(AD  461  744,  Litton  Sjrs- 
taaa,  HADES-12-65-1;  1965) 

■emlaphere  cone 

(rc/rb  :  0.07) 

0.3 

0.5 

0.65 

0.8 

0.9 

1.0 

1.1 

1.2 

1.3  .  . 

2  a  10* 

3.14  a  10* 

3.3  a  10* 

4.15  a  10* 

4.30  a  10* 

4.40  a  10* 

4.40  i :  10* 

4.38  a  10* 

_4.30  a  10* 

None 

Crahea,  T.  J.  and  Butlar, 

C.  B. 

(AD  240  084,  AEDC  TN  SO¬ 
US;  1960) 

Healaphere  cone 
cylinder  end  el- 
Upaolde  cylinder 

0.50 

0.  B0 

0.85 

0.90 

0.95 

0.975 

1.00 

1.05 

1.10 

1.20 

1.60 

Kangs  of  Is 

4.0  •  10* 
to 

7.2  •  10* 

None 

Eatz,  J.  I. 

(AD  159  801,  NAV0KD  15849; 

1958) 

Baalaphere  cylinder 

0.26 

0.49 

0.69 

0.83 

1.22 

1.6 

2.15 

2.48 

4.20 

(Rs  not  reported) 

Modified  Newtonian 

— - - I _ 

Hemisphere  cone 

3.12 

5.64 

5.80 
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TABLE  1.  (Continued) 


■  ■  —  —  — 

Sow <» 

HP151SS 

Ksch 

Nua6er 

ressramrasBrara 

IT*)' 

Theoretical 

Comparison 

Lord,  D.  1.  and  Ulnan.  E.  F. 

Oglva  Cylinder 

4.04 

Perturbation  the- 

(KACA  EH  L31L20;  1932) 

ory  end  method  of 

characteristics 

Im,  E.  f. 

BMlaphtrt 

1.97 

(Ba  not  raportad) 

Modified  Newtonian 

(Journal  of  Aaroapaea 

3.8 

I  Equation 

fclaocaa.  Hay  1*39) 

3.8 

6.9 

Cow  6  *  43* 

6.9 

Generalized  New- 

too lan  Equation 

Ogl  *.  [tf>/l  -  4) 

Exact  theoretical 

and  th(4/L)  •  0.3] 

methods  compared 

Com  0  -  <1" 

with  Generalized 

n 

Hew too lan  Theory 

Equation 

Hattkawa,  C.  U. 

Prolate  spheroid 

.«-■  Bjlj 

Frond tl-Clauret 

(MCA  n  2319,  193A 

S/d  •  6,10 

!  El 

correction  to  in- 

■rflS 

■  1 

compressible  po- 

tentlal  flow  solu- 

) 

0.92 

tlon 

0.94 

Ulchola,  J.  I. 

0.90 

Hone 

(AD  223  M2  M2.  AIDC- 

0.93 

4.10* 

j 

W-39-100) 

[•  ^  _ , 

1.00 

1.03 

} 

1.10 

1.20 

1.30 

1.40 

1.30 

| 

1.60 

Chaachl,  E. 

■auiapharlcally- 

(Ba  not  raportad) 

■onn 

(AD  291  399,  IMWI 

cap  pod  oglva 

!  KjLfl  . 

B  *048;  .1962) 

O/r^  -  0.16) 

(aaa  alao  Johanaoa,  V.  F., 

!i  A  « 

■OH  W  4061-771  1961) 

* 

1.33 

Van  Dyke 

2.17 

3.24 

Method  of  Charac- 

4.84 

terlstlcs 

Faaraon,  A.  0. 

Hemisphere  cone 

0.23 

3.0  a  10* 

Hone 

(MSA  n  D1961;  1963) 

0.40 

2.3 

0.60 

3.1 

0.60 

3. 

0.83 

3.6 

0.90 

3.9 

0.95 

4.0 

1.00 

4.1 

1.03 

4.1 

1.10 

4.2 

1.15 

4.2 

1.20 

4.2 

1.47 

3.0 

1.60 

3.0 

1.90 

3.0 

| 

Farklna,  S.  W. ,  Jorganaan, 

Heat sphere  cylinder 

rmnw 

Hew  ton  lan  Theory 

L.  ■.,  and  Somnar,  S.  C. 

Equation  1 

(BACA  E  1166;  1936 

ISJmJCB . 

Hemisphere  cone,  9  • 

1.50 

Cone  Theory 

6*39  *  c 

1.97 

3.10 

Hemisphere  cone,  0  - 

1.54 

rrmwanK^i; 

S’10'  e 

1.97 

3.06 

i; .  JKsgsr  jMSBXM' 

j 

Spraltar,  J.  B.  and  Alkana, 

Parabolic  arc  bodias 

0.80 

(la  not  raportad) 

Perturbation 

A.  T. 

S/d  >  6,  10,  and  6/2 

0.83 

theory 

(■ASA  TB-2;  1939) 

0.90 

1.0 

Cones 

1.0 

Stlna,  B.  A.  and  Uanlaaa, 

Baaiaphara  cylinder 

1.97 

3- 10*  to  6-10b 

Newtonian  Theory 

K. 

3.04 

2.6  •  10*  to 

Equation  1  and 

(CACA  T»  3344;  1934) 

4  •  10* 

Incompressible 

3.60 

2.64  •  10* 

Theory 

20 
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From  the  body  geometry,  L  may  be  determined  from  the  given  values 

and  R  . 
c 

sin  m  L/R 

*  - 
If 

r^  «  0.592  inch 

R  »  5.534  inches 
c 

then 

ip  m  27  deg 

and 

L  ^  2.50  inches 

0„  *  local  slope  at  the  nose  (x  »  0) 
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3.  Values  for  C  and  for  (P/P  )  may  be  found  In  Ref.  11,  or  in 
PN  '  C'N 

Che  National  Advisory  Committee  for  Aeronautics  publication,  NACA  Report 
1135,  Washington,  D.  C.,  1953;  for  example, 

C  -  0.53  and  /p/P  \  -  0.29 

'  Vn 


From  Eq.  6, 


C 


may. 


C  /sin20vt  -  2.59 
PN  N 


4.  Modified  Newtonian  theory  gives: 


x/D 

0 

Theory 

1.70 

C 

_ 2 _ 

Experiment 

1.68  (Stagnation  value.  Fig.  1) 

0.046 

0.46 

0.46 

0.325 

0.34 

0.36 

0.604 

0.24 

0.26 

0.883 

0.15 

0.16 

1.161 

0.084 

0.085 

1.440 

0.037 

0.030 

1.719 

0.012 

-0.03 

1.999 

0 

-0.08 

One  sees  that  the  modified  Newtonian  theory  begins  to  deviate  significantly 
from  experimental  data  as  6  becomes  small.  This  was  expected. 


5.  The  next  step  is  to  find  the  point  beyond  which  the  shock-expansion 
method  should  be  used.  Equation  10  may  be  used: 

(x'/L*)  -  0.60  +  0.15K 

m 

K  -  Kji/L)  -  2.17 
K  -  1.03 

(x'/LO  -  0.60  +  0.15(1.03)  -  0.75 
m 

x'  -  1.87 
m 
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x  -  1.74 
m 


(x/d)  =  1.47 

si 


6. 

theory : 


Compute  at  the  matching  point  using  modified  Newtonian 


C 

P 


sin2  9 


max 


mp 


C 

P 


*•» 


C  -  0.033 
P 


7.  From  the  definition  of  C  , 

P 


P  (y/2)PooMoo2 


Substituting 


P  ■  P  along  the  cone  surface 
w 


and 


I 


▼ 
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one  obtains 


p  \  -  >A  p. 

'■'•IT " 


(P/Pt) 

\  vN 


,cp +  p°°/q~K  +  *.a. 


From  these  equations,  the  shock-expansion  method  may  be  started. 

For  the  first  increment,  use  C  ■  0.033  at  x/d  ■  1.47  (modified  Newtonian 
value) . 


(P/Pt)1  ’  °'033  +  °'30  (  0.56+^0.30  ) 


rjico  -  2'y*2co  -  0-30 


From  a  table  of  Prandtl-Meyer ,  functions  for  P/Pt  -  0.112,  «  28.56. 

8.  Compute  the  Cp  at  the  next  station  downstream.  Station  2,  as 
follows : 


A9  ■  1 02  -  0 


02  -  vx  +  1 3 . 47  -  6.5 j  «  28.56  +  3.03  -  31.59 


WpSr"'SS  'jpi  * om 


From  the  tables. 


HWC  TP  4570 


(p/pt). 


0.095 


=  +p./9«, 


(p/pt)  TTi - P"/q" 

'  ,2  (p/pt)„ 


-  /«  /v«r\  0.56  +  0.30  n 

C  -  (0.095)  - pr^jo - 0.30 

p2 


C  -  0.28  -  0.30  -  -0.02 

P2 


Repeat  the  process  for  the  next  station.  Station  3: 


P/Pt  -  0.095 


v2  -  31.59 


v2  -  31.59  +  |e3  -  e2|  -  31.59  +  |0  -  3.47 1 


»  35.06 


K),  -  0. 


C  -  0.0763  (2.96)  -  0.30 


C  -  -0.074 
*3 
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9.  For  the  cylindrical  portion  of  the  body,  use  Eq.  lit 


p4 


C  exp(-Ax/LK) 


Ax  Ax/p  2.337  -  1.999 
L  "  L/p  "  2.111 

K  “  1.03 
Ax  0.338 


LK  2.111(1.03) 
exp(0.155)  -  1.167 
„  _  -0.074 


0.155 


p4 


1.167 


-0.063 


C  -  -0.074  exp (-0.658/2. 175)  -  yr-Py4 

pr  l.JjJ 


C  *  -0.054 
P5 


f  this 

x/D 

process  nay  be  summarized 

C 

P 

as  follows 

Calculation 

Experiment 

0 

1.70 

1.68 

0.046 

0.46 

0.46 

0.325 

0.34 

0.36 

0.604 

0.24 

0.26 

0.883 

0.15 

0.16 

1.161 

0.084 

0.085 

1.440 

0.037 

0.030 

1.47 

0.033 

1.719 

-0.020 

-0.025 

1.999 

-0.074 

-0.060 

2.337 

-0.063 

-0.050 

2.657 

-0.054 

-0.040 

The  calculations  and  data  appear  in  Fig.  7. 
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EXAMPLE  2:  EQUATION  15 

As  an  example  of  the  use  of  £q.  15,  consider  an  ellptlcally  nosed 
cylinder  (see  Fig.  13)  R/L  ■  0.61  at  II  *  0.90i 


Cp  -  cp  sin28b  +  tt/«F(eb,  *.) 


max 


1.  From  Fig.  1,  C  *  1.215  at  M  *  0.9. 


max 


2.  For  these  coordinates. 


tan  6  *  (R2/L2) (-x/r) 


where 


r  -  %/(R2/L2  )  (L2-x2) 


-X 

b,  deg 

p 

exp 

2.684 

68 

1.075 

2.359 

45 

0.42 

2.033 

33.5 

0.05 

1.464 

20.7 

-0.15 

0.935 

12.3 

-0.58 

0.488 

6.25 

-0.75 
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3.  For  convenience  In  computation,  choose  calculated  points  at 


9b 

sin20^ 

C  sln20, 

P  b 

max 

F(V  M.) 

(R/L)F^0b, 

90 

1.000 

1.215 

0.0 

0 

80 

0.9698 

1.18 

-0.037 

-0.024 

70 

0.8830 

1.072 

-0.088 

-0.045 

60 

0.7500 

0.911 

-0.165 

-0.10 

50 

0.5868 

0.714 

-0.276 

-0.168 

40 

0 .  ( 131 

0.502 

-0.422 

-0.257 

30 

0.2500 

0.304 

-0.600 

-0.368 

20 

0.1170 

0.142 

-0.798 

-0.486 

10 

0.0301 

0.037 

-1.00 

-0.61 

0 

0 

0 

-1.187 

-0.724 

eb 

C 

P 

Calculation 

Experiment 

90 

1.215 

80 

1.133 

70 

1.00 

68 

1.075 

60 

0.800 

50 

0.536 

•  •  •  •  • 

45 

0.42 

40 

0.235 

33.5 

0.05 

30 

-0.068 

20.7 

-0.15 

20 

-0.345 

12.3 

-0.58 

10 

-0.574 

6.25 

-0.75 

0 

-0 . 724 

These  results  are  presented  in  Fig.  13 
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PERTURBATION  THEORY 


For  some  bodies  (pointed  bodies  or  those  which  are  treated  in  terms 
of  equivalent  pointed  bodies,  in  the  transonic  regime  for  instance)  and 
for  Mach  numbers  other  than  those  specified  above  as  being  ammenable  to 
the  modified  Newtonian  or  the  modified  Newtonian/ shock-expansion  methods 
(M^  <  0.7),  the  most  useful  methods  for  solution  are  those  of  small 
perturbation  theory  as  developed  by  Van  Dyke  (Ref.  21,  22,  and  23)  and 
Sprelter  (Ref.  24  and  25).  Somewhat  simpler  linearized  theories  have 
been  developed  by  Laitone  (Ref.  26  and  27)  and  Adams  and  Sears  (Ref.  28). 


Several  texts  (Ref.  29  and  30)  contain  explanations  of  and  background 
material  on  perturbation  theory.  Van  Dyke  (Ref.  21  and  22)  presents  more 
detailed  developments.  The  pertinent  developments  of  the  theory  are  fully 
presented  in  the  cited  references;  only  the  ?quations  directly  employed 
in  the  computation  of  the  pressure  distribution  will  be  presented  here. 


The  application  of  perturbation  theory  begins  with  the  description 
of  the  flow  over  a  body  by  means  of  a  velocity  potential.  This  approach 
presumes  irrotatlonal  flow.  This  idealization,  however,  seems  to  be 
justified  by  the  results  in  the  subsonic,  transonic,  and  low-supersonic 
regions  obtained  by  Van  Dyke,  Sprelter,  and  others  (see  Fig.  8,  14,  and 
15,  and  Ref.  21-30). 


FIG.  15.  Pressure  Distribution  on  a 
Parabolic  Arc  Body,  £/d  -  6. 
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In  this  method,  the  basic  flow  field  is  described  by 

U,  -  u 
1  00 


u2  -  0 


u3  *  0 


and  when  the  presence  of  a  body  disturbs  the  flow  field,  by 


U,  “  u  +  u. 
1  00  1 


U2  "  u2 


U3  “  u3 


are  the  perturbation  velocities.  The  problem  Is 
of  the  nondimens lonallzed  parameters, 

Uj/Uoo*  u2/uoo’  u3/u« 

The  criterion  for  application  of  this  method  is  that  these  nondimensional 
parameters  are  much  less  than  unity.  For  the  majority  of  the  bodies  of 
interest,  this  criteria  will  be  met  except  In  the  region  of  the  stagnation 
point  of  blunt  bodies.  Van  Dyke  (Ref.  21  and  23)  presents  a  method  of 
correction  for  subsonic  flow  that  renders  the  solution  valid  even  in  the 
stagnation  region  if  the  nose  contour  is  no  more  blunt  than  round.  A 
rather  obvious  example  of  a  body  that  would  be  excluded  from  treatment 
by  perturbation  methods  is  a  flat- faced  cylinder.  Such  a  body  is  excluded 
because  its  nose  is  more  blunt  than  round,  and  the  ratios  u^u^  and  u-j/u^ 
would  not  be  small. 

The  first  step  in  the  development  of  perturbation  methods  is  the 
description  of  the  flow  field.  Both  velocity  potential  and  stream 
functions  have  been  used  for  this  purpose.  /.  perturbation  velocity 
potential,  <f>,  is  most  often  used  and  defined  such  that 


where  u_,  u2,  and  u. 
then  stated  in  terms 


....  ‘ 

.  r  '  1  V- 


-*?<••**  if# iVWlSWWW 


3x  yx  u 


3y  y  *»« 


M-*  -S 


9z  Yz  u 


The  equation  of  motion  of  the  flow  field  about  the  body  may  then  be 
stated  in  terms  of  this  perturbation  velocity  potential: 


( 1  —  M  2 )  4>  +<js  +<j) 

°°  XX  Yyy  Yzz 


M  2  [  /  2<j>  +<J)2  +  ({)2  +  (j)  2\/$  +0  +  <j>  \ 

°°  l  2  \  Yx  Yx  Yy  Yz  /\  xx  yy  Yzz J 


+  2(j)  $  +  4>  24>  +  ct>  24>  +  <J>  2(J)  +  2<t>  /l  +  <j>  U 

x  xx  x  xx  y  yy  z  zz  z  y  x 12 


+  2*y(l  +  *x)+xy] 


Equation  16  will  yield  velocity  components  which  are  exact,  as  long  as 
the  limitation  that  <Px,  <}>y,  4>z  <<:  1  is  met.  These  equations,  however, 
are  quite  unmanageable.  Converting  to  a  cylindrical  coordinate  system, 
reducing  the  equations  to  second  order  in  terms  of  the  derivatives  of  <J>, 
and  limiting  consideration  to  bodies  of  revolution  yields  the  following 
form: 


<p  +  —  ±  3  24>  =  m  2  2 - ±  1  3Z4>  <P 

rr  r  xx  °°  2g2  Tx  xx 


(y+l)M 


+  2d>  <t>  +  <bz  <b 

r  xr  r  rr 


for  subsonic 


for  supersonic^ 
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where 


6 


subsonic  flow 


supersonic  flow 


A  linearized  form  for  determining  the  potential,  <f>,  is  obtained  by 
further  approximation,  neglecting  the  terms  on  the  right-hand  side  of 
Eq.  17  which  contains  products  of  the  derivatives  of  <J>.  This  then 
results  in  the  linear  equation 


4>  + 

Trr 


—  ±  -  0 
r  xx 


This  linearized  form  has  been  used  by  Laitone  (Ref.  26  and  27)  and  by 
Adams  and  Sears  (Ref.  28)  in  both  the  subsonic  and  supersonic  flight 
regimes.  Initially  the  coefficient  of  pressure  resulting  from  linearized 
theory  was  expressed 


Cp  -  -2*x  (19a) 


However,  as  pointed  out  in  Ref,  28,  <J>X  and  <PT  are  actually  of  different 
order  at  the  body  surface,  and  the  accuracy  of  the  method  in  predicting 
the  coefficient  of  pressure  is  improved  by  retaining  the  second-order 
term  ir.  <j>r, 

C  -  -2<p  -  <p  2  (19b) 

p  x  r 

This  improvement  comes  at  no  increase  in  complexity  of  the  solution,  for 
at  the  body  surface 

-  dR/dx  (20) 

where 

R  *  R(x)  describes  the  variation  of  the  cross-sectional 
radius,  R,  with  axial  distance,  x. 
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It  should  be  noted  that  Matthews  In  Ref.  31  employed  Eq.  19a  rather  than 
Eq.  19b  for  linearized  theory.  For  this  reason  his  statement  in  that 
reference  as  to  the  Inaccuracy  of  linearized  theory  is  not  correct  as 
is  pointed  out  by  Adams  and  Sears  in  Ref.  28.  Adams  and  Sears  (Ref.  28) 
presented  the  following  approximate  expression  for  the  linearized 
potential  in  the  subsonic  flow  regime: 


4>(x,r)  «  ~  [s’  (x) In  \  S'(0)ln  x  -  \  S'(Jt)lnU-x)  (21a) 


-  \  f*  S"(C)ln(x-C)dC  +  S"(5)ln(C-x)d? 


where 


S  -  ttR2 
R  -  r(x) 

l  -  total  length  of  body 


Notice  that  this  solution  yields  singularity  at  the  stagnation  point; 
for  this  reason,  the  solution  obtained  is  not  valid  In  the  region  of  a 
stagnation  point.  However,  for  points  a  distance  x  >  0.05£  away  from 
a  stagnation  point,  this  solution  gives  good  results.  For  pointed 
bodies,  the  solution  is  valid  except  very  near  the  tip.  Additionally, 
because  Eq.  21a  was  developed  for  subsonic  flow,  this  solution  presumes 
that  the  body  curvature  is  smooth  and  continuous  f-rom  nose  to  tall. 

This  would  eliminate  from  consideration  those  configurations  with 
cylindrical  afterbodies  or  flat  stems.  However,  comparison  of  the 
results  of  Eq.  21a  with  experimental  data  for  configurations  with 
cylindrical  afterbodies  shows  that  the  solution  is  good  except  in  the 
region  of  discontinuity  (Ref.  24  and  25).  The  corresponding  supersonic 
flow  solutions  developed  from  linearized  theory  are  limited  to  pointed 
bodies  but  the  shape  of  the  aft  section  has  no  influence.  As  a  general 
rule,  approximate  equations  such  as  Eq.  21a,  developed  from  the  linearized 
equation,  Eq.  18,  will  obtain  good  results  for  <  0.7.  As  an  example, 
the  pressure  distribution  over  a  prolate  spheroid  Tan  ellipsoid  of 
revolution)  in  subsonic  flow  is  calculated  in  the  Example  of  Linearized 
Perturbation  Theory,  page  37.  Equation  21a  yields  the  coefficient  of 
pressure  at  the  body  surface: 


C 

P 


4t2 

i2 


(21b) 
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where 


R  -  2t II  (ix  -  x2)*5 
l  ■  total  length  of  body 
T  “  d/i 

The  Cp  at  the  stagnation  point  is  again  given  by  Eq.  13. 

Solutions  obtained  from  the  linearized  equation,  Eq.  18,  cannot  be 
used  in  the  transonic  regime  because  the  term  82$xx  retained  is  no  longer 
large  in  comparison  to  one  of  the  neglected  terms  that  appeared  on  the 
right-hand  side  of  Eq.  17.  To  correct  this  difficulty.  Van  Dyke  (Ref.  21) 
and  Spreiter  (Ref.  24  and  25)  present  more  rigorous  approaches  to  the 
solution  for  the  perturbation  velocity  potential  that  offer  considerably 
improved  accuracy,  particularly  in  the  transonic  flight  regime.  This 
improvement  is  obtained  by  including  a  nonlinear  term  on  the  right-hand 
side  of  the  potential  equation,  Eq.  17: 


(1  -  Kj)  ^  +  (1/r)  <J>r  +  <|>rr  -  Ha>2(Y+l)^xx  (22) 


This  equation,  for  a  frictionless,  lrrotatlonal  flow,  provides  an  improved 
first-order  solution  which  is  applicable  to  the  subsonic,  transonic,  and 
supersonic  regimes.  Again,  the  equation  for  determining  the  coefficient 
of  pressure  at  the  body  surface  is 


C 

P 


Van  Dyke  (Ref.  21)  presents  a  solution  to  Eq.  22  and  develops  a 
correction  method  that  renders  the  solution  valid  in  the  vicinity  of  the 
stagnation  point  of  blunt  bodies  in  subsonic  flow,  with  the  limitation  on 
the  solution  being  only  that  the  ends  be  no  more  blunt  than  round.  The 
basic  equations  of  Van  Dyke’s  solution  may  be  presented  as  follows: 

MC02(Y+1) 

4>(x,r)  -  (F  +  f)  ln(r)  +  (G  +  g)  + - F(F'  ln(r)  +  G’)  (23) 

2(Mro2  -  1) 
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where  the  prime,  * ,  denotes  differentiation  with  respect  to  x,  and 


F(x)  -  RR' 


f(x) 


G(x)  =  { 


(Y+l) 

~  1) 


F(x)  In 


FF’ 


8 


2/(x-a)  (b-x) 


+  2 


1  F(x)  -  F(£)  ,r 

?Ja  l*-^l  ? 


.  '  1 


FW  ln  Toby + f*  ae  ,  m„<i 


g(x)  - 


f(x)  ln 


8 


2/(x-a)  (b-x) 


X  fb  f(x)  -  f(Q  M 

2 ja  dt  ■  H« 


<  l 


f(x)  ln 


8 


2(x-a) 


/. 


ilzLz  UP.  dr 

X  -  £  ^ 


,  M  >  1 

9  oo 


As  a  correction  factor  for  blunt  ends  in  subsonic  flow.  Van  Dyke  then 
develops  a  correction  factor  to  be  added  to  Eq.  23: 

C  .[  1  F*(a)  1  F2(b)  1  /(X-X)  \ 

R  I  2  x-a  2  b-x  I  \4(1  „ jJ 


which  renders  the  solution  valid  in  the  region  of  the  stagnation  point. 
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As  an  example  of  this  method,  Van  Dyke  conjiders  the  case  of  a 
prolate  spheroid  (Ref.  21)  and  arrives  at  the  result  for  C  at  the 
body  surface:  P 


C 

P 


1  -  5x2  +  2x4 
(1-x2)2 


Adding  the 


correction  factor  for 


the  round  ends  gives  a  final  expression 


C 

P 


2 


+  1  + 


/  Y+l  \ 

1  -  llx2  +  4x4 

1  4  '  (1  -  M  2) 

OO 

(1-x2)2 

Van  Dyke  used  the  following  illustrated  coordinate  system  for  his 
development: 


Results  obtained  from  Van  Dyke's  equation  are  compared  with  those  of 
linearized  theory  in  the  example  presented  on  page  37,  for  a  prolate 
spheroid  T  *  1/6  at  *  0.70.  This  comparison  indicates  that  for 
this  Mach  number  there  is  no  advantage  in  using  Van  Dyke's  equation, 
for  the  improvement  which  results  is  very  little,  if  any,  except  in 
the  stagnation  region.  The  use  of  Van  Dyke's  method  is  not  so  straight¬ 
forward  as  is  linearized  theory,  Eq.  21a.  However,  in  the  transonic 
regime  the  accuracy  of  Van  Dyke's  method  improves,  whereas  that  of 
linearized  theory  degenerates.  The  results  from  Van  Dyke's  second- 
order  method  are  given  in  Fig.  14  for  ■  0.90  and  x  ■  1/6. 

Spreiter  and  Alksne  (Ref.  24)  present  a  very  complete  development 
for  the  solution  to  Eq.  22  for  pointed  bodies  of  revolution.  These 
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results  are  illustrated  by  Fig.  15,  which  has  been  reproduced  from 
Ref.  24.  Again  the  process  of  solution  is  rather  difficult  but  for 
pointed  bodies  in  the  transonic  regime,  there  is  little  alternative. 

For  the  details  oi  this  solution,  the  reader  is  referred  to  Ref.  24. 

Van  Dyke  (Ref.  21)  also  presents  a  complete  second-order  solution. 
In  this  solution,  all  of  the  terms  on  the  right-hand  side  of  Eq.  17  are 
retained  but  are  handled  as  follows: 


0r  ,  ,  f(Y+D  ] 

0„  +  ±  8Z0 _  =  M  2 (2  -  ±  I  B2 


rr  r  xx  00 


—  1  1  Mx* 


+  2$  <f>  +  4>  2<j> 

rxr  r  rr| 


for  subsonic  \ 
for  supersonic^ 


where 


1  -  M_ 


subsonic  flow 


supersonic  flow 


The  complexities  of  such  a  solution,  however,  carry  it  beyond  the  scope 
of  this  report,  and  mention  of  this  method  is  made  only  for  the  sake  of 
completeness.  Okauchi  (Ref.  12)  compares  the  results  of  Van  Dyke's 
second-order  solution  for  an  ogive  with  data  obtained  from  tests  of  the 
hemispherically  capped  ogive  (rc/rb  ~  0.15),  discussed  previously. 
Representative  results  for  low  supersonic  Mach  numbers  are  included  in 
Fig.  8. 


EXAMPLE  OF  LINEARIZED  PERTURBATION  THEORY 

As  an  example  of  perturbation  methods,  a  prolate  spheroid  with 
T  =  1/6  at  Mot  *  0.70  is  considered.  At  the  surface  of  a  prolate  spheroid, 
the  dimensionless  radius  is 


R(x)  =  (2x /l)  (xSL  -  x2)*5 


T  =  d/SL 


37 


NWC  TP  4570 


At  the  surface,  r  -  R(x);  the  cross-sectional  area  is 


S(x)  -  R2  (x)  -  —  (xi,  -  x2) 


S'  (x)  -  (£  -  2x) 


S"(x)  -  - 


8ttt' 


From  Eq.  21a 


In  (*§)  -  \  S’  (0)  In  x 

-  f  S'  (£)  ln(2  -  x)  -  S"(£)  ln(x  - 

+  k  fl  S"CO  ln(C  -  x)ds| 
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for  small  values  of  r  (i.e.,  in  the  vicinity  of  the  surface).  Mike  the 
substitution  for  S’  and  S": 


<Kx,r)  *  ^  a  -  2x)  In  (^|)  -  | 


In  x  +  —  In  (£  -  x) 


+  /  ln(x  -  £)d£  -  /  ln(£  -  x)d£l 


Differentiate  with  respect  to  x: 


2t2 

4>v(x,r)  =  - 

x  a2 


(*)-!(&) 


^ j  dC  +  ln(0)  -  o) 


d£  +  0  -  ln(0) 


~  2t2 

K  -  ■—  -2  m 

x  s,2 


^  -  srbr  •/“*" 


■r* 


+  I  d£  +  2  ln(0) 


4(ix-x2) 


+  -j  (ln(x-£)}’ 


-  |[ln(5-x)]*  -  In  (0) 
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p  z  — 
x  i2 


m  (if) 


jz  1 

+ - -  —  ln(£x-x2 ) 

4(i,x-x2) 


Now  evaluate  <p  at  the  body  surface  r  -  R(x)  -  (2t/H)  (Jlx-x2)*5: 


-4t2 

<J>  (x,R)  s  -  In 

X  Z2 


(1) 


+  (2t/Z)  (Zx-x2) 

4(S,x-x2)  (2t/Z) 


Combine  the  first  and  last  terms  using  the  definition  of  R: 


=  =*£  In  (If)  +  — 

X  £2  [  21  4(£x-x2) 


Also  at  the  body  surface, 


Ux-x2)^ 


From  Eq.  19,  the  coefficient  of  pressure  is 


C  =  -2<J>  -  <p  2 

p  x  r 


p  i2 


p  n2 


(H) 


(t2/£2) (A-2x) 2 


4(Jlx-x2)  (£,x-x2) 


(^f) 


i  r  -  4xJi 


4(£x-x2)  L  £x-x2 


+  4x2 ]  I 

— li 


{•  ■■  R) 


+  1  + 


4(Jlx-x2) 
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The  last  terra  of  this  equation  in  terms  of  the  body  cross-sectional 
radius  is  expressed 


C 

P 


2  In 


+  1  +  — 


This  may  be  compared  with  Van  Dyke's  solution  (page  42)  by  substitution 

of  £  *  2  and  R  =  Tv/l  -  x2  to  fit  Van  Dyke's  choice  of  coordinates.  Thus 
the  linear  theory  gives 


C 

P 


In  +  1  + 


1  -  x' 


as  compared  to  Van  Dyke's  solution 


C 

P 


2  xn  —  +  1  + - — - 

1  -  x‘ 


Y+l  T4  1  -  llx2  +  4x4 
1  -  M^2  (1  -  x2)2 


The  difference  between  the  two  equations  is  in  the  second  bracket. 
This  arises  from  Van  Dyke's  use  of  Eq.  22  and  the  correction  term 
instead  of  3q.  which  is  the  starting  point  for  the  development  of 
Eq.  21a. 

As  a  numerical  example,  consider  a  prolate  spheroid  with  x  ■  1/6 
at  ■  0.7  (data  from  Ref.  31). 
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From  these  calculations,  one  sees  that  the  more  complicated  methods 
of  Van  Dyke  offer  little  advantage  at  -  0.70.  However,  it  should  be 
noted  that  at  Mach  numbers  in  the  transonic  regime,  the  accuracy  of 
linearized  solution  degenerates,  whereas,  the  results  obtained  by  Van  Dyke* 
method  improve  in  accuracy  somewhat  (see  Fig.  14). 


SIMILARITY  RULES 

The  final  method  to  be  discussed  is  that  of  similarity.  In  this 
method,  a  known  pressure  distribution  on  one  body  is  used  to  predict 
the  pressure  distribution  on  another  affinely  related  body.  Affine 
Implies  that  all  coordinates  in  a  given  direction  change  only  by  a 
uniform  ratio  from  one  member  to  another.  For  example,  on  bodies  of 
revolution  this  affine  relationship  may  be  expressed 

r  -  tR(x) 

where  r  is  the  radial  distance  to  the  surface  from  the  axis  of  symmetry, 
and  T  is  a  characteristic  body  parameter  such  as  fineness  ratio  or  slope. 

This  approach  may  be  used  in  two  ways: 

1.  In  subsonic,  transonic,  supersonic,  oi  rsonic  flow,  one  may 
use  experimental  measurements  on  one  member  of  an  affinely  related  family 
of  bodies  to  predict  the  distribution  on  the  other  members  of  the  family. 

2.  In  subsonic  flow,  one  may  correct  exact  incompress  id  potential 
flow  solutions  for  the  effects  of  compressibility. 

Similarity  rules  are  discussed  and  compared  with  data  in  Ref.  .20,  29,  30, 
31,  and  32. 

For  irrotational  flows,  a  velocity  potential  exists  such  that 

$  -  <t>(x,r;  M,y,t) 

The  solution  (Ref.  33)  proceeds  by  separating  the  dependent  variable 
into  components  and  then  stretching  each  component  and  the  independent 
variables  by  factors  that  depend  on  the  original  parameters: 

4/u^  -  x  +  0  +  $ 
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where  4>  is  the  linearized  potential  theory  and  0  the  second-order 
increment.  For  the  linearized  problem 


(1  -  M 


!)<t>  + 

xx 


rr 


+  0 


(25) 


<p  -*■  0  at  infinity 


4>r  -*■  xR'(x)  at  r  ■  xr(x) 


Thus,  the  linearized  problem  is  not  dependent  on  y  and  depends  only  on 
M  and  T.  The  similarity  rules  are  obtained  by  considering  arbitrary 
scale  transformation  in  <p  and  r.  The  only  choice  that  reduces  the  number 
of  parameters  in  the  problem  is 

<J>(x,  R)  -  1/32  F(x,  p)  (26) 


P  -  3r 


B 


for  subsonic 


for  supersonic 


The  problem  then  becomes 


F  +  F  <{)  ±  F  =0  (27) 

^  XX  /+  for  subsonic  \ 

\  -  for  supersonic/ 


Thus,  the  original  problem  is  reduced  to  one  where  dependence  is  on  Mro 
and  x,  not  separately,  but  only  in  the  combination  of  8x.  The  resultant 
similarity  rule  is  then 


4>(x,r;  M^.x)  =*  (1/B2)  F(x,6r;$x) 


(28) 
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For  the  velocity  potential  and  the  coefficient  of  pressure  to  its  first 
order,  the  results  are 


(l/uoo)4>(x,r;  M^.t)  -  x  +  (l/B2)F(x,8r;8x) 


(29) 


C  (x,r;  M,t)  -  (l/82)P(x,8r;BT) 
P 


(30) 


Two  further  equations  used  for  determining  the  Cp  in  compressible 
flow  from  the  incompressible  C  are  stated  by  Sears  (Ref.  20): 


inc 


C  (x.R)  -  C  (x,R)  -  (R2)"  In  6 
P  pinc 


(31) 


and  by  Lees  (Ref .  34)  : 


C 


P 


C 

Pinc 


In  g 

0. 31(ln  t) 


Sears  (Ref.  20)  indicates  that  for  slender  bodies  of  revolution  in  low 
subsonic  flow,  <  0.7,  correction  formulas  such  as  these  closely 
approximate  theoretical  results.  This  conclusion  was  reached  by  con¬ 
sideration  of  the  results  from  Ref.  31  for  ellipsoids  of  revolution 
with  X,/d  ratios  of  6  and  10. 

In  the  transonic  flow  regime,  the  similarity  rule  takes  a  somewhat 
different  form.  Additionally,  for  bodies  of  revolution  in  the  neighbor¬ 
hood  of  the  body  surface,  a  similarity  rule  exists  only  if  the  body  is 
very  slender.  In  terms  of  determining  the  pressure  distribution,  this 
severe  limitation  makes  the  use  of  transonic  similarity  rules  of  little 
value  (Ref.  32).  However,  if  the  bodies  of  interest  are  sufficiently 
slender,  affinely  related  bodies  will  be  similar  if  the  similarity 
parameter, 

J  «  T2  In  [/2  t(1  -  Mj]  (32) 

is  equal  for  both  bodies. 
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COMPUTER  TECHNIQUES 

Although  consideration  of  complex  computer  techniques  for  determining 
the  pressure  distribution  was  not  the  objective  of  the  program,  It  was 
noted  during  the  literature  survey  that  several  computer  techniques  are 
being  used  with  good  results  for  blunt  bodies  In  supersonic  flow;  some 
of  these  are  noted.  A  brief  comparison  of  theoretical  and  experimental 
results  is  given  in  Ref.  35.  Methods  developed  by  NASA  (Ref.  36)  and 
at  NWL  (Ref.  37)  are  representative  of  those  that  yield  good  results. 

Also,  Ref.  38  reports  a  modification  of  a  North rop-developed  program  that 
yields  good  results  in  the  Mach  1  to  3  range. 
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Section  3 


HEAT-TRANSFER  DISTRIBUTION  ON  BLUNT  BODIES 
OF  REVOLUTION 


The  analysis  of  the  convective  heating  of  a  high-speed  vehicle  is 
an  Integral  part  of  the  evaluation  of  the  aerodynamic  and  thermal 
behavior  of  the  vehicle.  Methods  of  analysis  for  convective  heating 
have  been  extensively  developed  over  the  past  two  decades.  However, 
this  work  has  been  directed  primarily  toward  an  understanding  of  heating 
conditions  where  the  stagnation  enthalpy  is  much  larger  than  the  enthalpy 
of  the  air  at  the  surface  of  the  vehicle. 

The  experimental  data  reviewed  were  those  obtained  for  Mach  numbers 
less  than  5  both  in  free  flight  and  in  wind  tunnel  tests.  Not  all  of 
the  available  data  are  reviewed  since  some  wind  tunnel  conditions  were 
not  compatible  with  actual  flight  conditions  expected  at  the  corresponding 
Mach  numbers,  and  some  data  were  originally  presented  in  such  a  fashion 
that  it  made  comparison  difficult  to  achieve. 

Reasonable  agreement  between  theory  and  experiment  appears  to  have 
been  achieved  for  flight  regimes  above  a  Mach  number  of  1.6  for  blunt 
body  vehicles.  However,  adequate  experimental  verification  of  the 
analytical  methods  developed  for  the  high  Mach  number  regimes  has  not 
been  obtained  for  the  range  of  flight  Mach  numbers  from  0.7  to  1.6. 

Since  the  Reynolds  number  encountered  in  the  low  altitude  (sea  level 
to  50,000  feet),  the  high-speed  flight  regime  under  consideration  is 
expected  to  be  high;  a  large  portion  of  the  boundary  layer  is  expected 
to  be  turbulent.  It  will  be  necessary,  therefore,  to  examine  methods  . 
of  analysis  for  the  laminar,  transitional,  and  turbulent  heat- transfer 
regimes  and  to  compare  these  methods  with  the  appropriate  experimental 
results. 


STAGNATION-POINT  HEAT  TRANSFER 

Various  methods  of  computing  the  stagnation-point  heat  transfer  to 
blunt  body  supersonic  vehicles  have  been  presented  in  the  literature. 
Probably  the  best  known  of  these  methods  are  those  of  Lees  (Ref.  39), 
Fay  and  Riddell  (Ref.  40),  Detra,  Kemp,  and  Riddell  (Ref.  41),  Vaglio- 
Laurln  (Ref.  42),  and  Cohen  (Ref.  43).  These  methods  were  developed 
to  predict  the  heat  transfer  to  blunt-body  vehicles  as  vehicle  speeds 
were  pushed  to  higher  and  higher  values.  A  major  objective  of  these 
studies  has  been  to  account  for  the  thermochemistry  effects  that  arise 
as  a  result  of  the  higher  stagnation  temperatures  encountered. 


'  ity  aw  A  ft:  i  cffTi 
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The  equation  derived  by  Sibulkln  (Ref.  44)  for  the  heat  transfer  at 
the  stagnation  point  of  a  sphere  has  been  used  by  many  investigators  to 
evaluate  the  heat  flux  at  the  stagnation  point.  Starting  with  Sibulkin's 
equation 


*5  „  0.4 


Ni^  -  0.763L  (B)^  (p e/\iey  Pr 


where  Nu^  is  the  Nusselt  number  based  on  a  typical  length,  L,  and  B  is 
the  velocity  gradient  at  the  stagnation  point,  the  heat  transfer  at  the 
stagnation  point  may  L  written 


\  -  °'763  (^e)’5  (Pr)°'4  k(\  '  T„) 


This  expression  may  be  modified  to  read 


d  -  0.763  (Pr 


’■6)  (w)*  (\  ■  h“)  <b)! 


Substituting  in  the  expression  for  the  velocity  gradient  at  the  stagnation 
point. 


£). 


one  obtains  the  result 


0.763  Pr 


'°'6)  (w)”  (\  '  h») 


where  Pr  ■  0.71  fcr  computations.  In  this  expression  the  density,  p  , 

s 


and  the  viscosity,  y  »  are  evaluated  at  the  outer  edge  of  the  boundary 

e 

s 


layer  at  the  stagnation  point  of  the  body.  Note  that  the  only  effect  of 
the  wall  temperature  in  this  equation  occurs  in  the  evaluation  of  the 
wall  enthalpy. 
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To  overcome  this  deficiency,  Fay  and  Riddell  (Ref.  40)  obtained  a 
correlation  of  numerical  results  with  the  following  expression  for  the 
heat  transfer  at  the  stagnation  point  of  a  blunt  body 


0.763 


H  (*->>'  (v-}J 


(38) 


The  velocity  gradient  at  the  stagnation  point  may  be  evaluated  from  the 
Newtonian  pressure  distribution: 


(39) 


Hence,  using  the  perfect  gas  equation,  a  heat-transfer  parameter  may  be 
written  as 


e 


Cohen  (Ref.  43)  has  also  presented  a  correlation  of  the  results  of 
numerical  solutions  for  high-velocity  heat  transfer  (U^  >  10,000  fpj) 
in  the  form 
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In  addition  to  these  expires s ions ,  Detra,  Kemp,  and  Riddell  (Ref.  41) 
have  presented  a  correlation  for  high  velocity  (U^  >  10,000  fps)  in 
the  simplified  form  of 


The  experimental  blunt-body  heat-transfer  data  available  for  com¬ 
parison  with  the  preceeding  methods  of  calculation  are  limited  for  flight 
Mach  numbers  between  2  and  5,  whereas  appropriate  data  below  a  Mach  number 
of  1.6  are  simply  not  available  (Ref.  45,  46,  47,  48,  and  49).  This  is 
shown  in  Table  2,  which  presents  a  summary  of  the  appropriate  blunt-body 
heat-transfer  distribution  data  available.  The  Mach  number  and  Reynolds 
number  at  which  the  data  were  obtained  are  indicated  in  the  table. 
Appropriate  stagnation-point  heat-transfer  data  were  obtained  from  the 
results  presented  by  Buglla  (Ref.  50),  Chauvin  (Ref.  51),  and  Chauvin 
and  Maloney  (Ref.  52).  Reference  should  be  made  to  the  Annotated  Bibli¬ 
ography,  page  87  in  the  Appendix,  for  an  indication  of  the  difficulties 
encountered  in  using  data  presented  in  the  other  reports  listed. 


<w 


\/  R_/P_  - 


2.28  x  10 


-2 


f 


TABLE  2.  Low  Supersonic  Blunt-Body  Heat-Transfer  Distribution  Data 
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m 

V 
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■ 

■  ■ 

m 
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m 

Local 
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2-5 

18.20 
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2 

X 

night 

Cappad  Cona 

2.47 

7.66  x  10* 

X 

X 

X 

J.  J.  lugli a 

In/**  *  0.7* 

2.0 

7.84  x  10s 

X 

X 

X 

Cona  half- 

2. B0 

8.36  x  10‘ 

X 

X 

X 

1 

angle  -  1*  J* 

2.97 

8.54  x  J0‘ 

X 

X 

X 

tg  •  8.491" 

3.14 

Naval  Ordnanca  Lab- 

Hind 

Healepherlcel 

2.80 

Not  reported 

0.366 
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X 

X 

X 

oratory 

Tunnal 

I.  Korobkln 

D  ■  1  In. 

Navord  Report  2863 

Naval  Ordnanca  Labo- 

Wind 

Haadaphara 

V  1  T^Bi 

Not  repotted 

0.866 

2.6  x  10J 

X 
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tunnal 

cylinder 

1.4  x  10s 

2 

K 

4.0  x  105 

X 

I.  Korobkln 

D  •  2  In. 

4.20 

4.1  x  105 

X 

Navord  Kapur t  3841 

4.87 

2.4  x  10s 

X 

Aaao  Kasaarch  Cantar 
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1.97 

20 
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2.4  x  10s 

X 

X 
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Haalaphara 

3.414 

n 
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HI 

ij 
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EB 
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Wa 

Ola 

mi  '  * 

mm 
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Wind 
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3.12 

16 

2.94 

* 

X 
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j 

N.  8.  DlaconiSg  K.  J* 
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3.12 
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X 

X 

Wlanlavakl,  and  J.  K. 
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D  •  1.4  in.  1 

NACA  -TO-4099 
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2.00 
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3.5 

X 

X 

X 

tart  Z.  E.  Backvlthg 

Tunnal 

D  -  3*|  In. 

!  IWtri  v  ■*. 

J.  J.  Callaghar 

Htai  >  ph*  rt- cy  1  *  .ida 

r 

1.57 

nQMi 

X 

X 

X 
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It  would  appear,  at  this  point,  that  the  results  of  the  experimental 
program  reported  by  Buglia  (Ref.  50)  represent  the  most  complete  set  of 
experimental  data  obtained  for  free-f light  or  wind  tunnel  environments. 
For  this  reason,  the  set  of  results  obtained  by  Buglia  has  been  employed 
for  comparison  with  the  formulas  previously  presented. 


In  the  flight  test  conducted  by  Buglia,  heat-transfer  and  transition 
data  on  a  highly  polished  hemisphere  cone  were  measured  over  a  Mach  number 
range  froth  2.32  to  3.14  up  to  a  free-stream  Reynolds  number  of  24  x  10^ 
based  on  the  body  diameter.  The  hemisphere  radius  was  6.498  inches  and 
the  cone  half-angle  was  14.5  degrees.  The  stagnation-point  heat-transfer 
results  obtained  from  this  experiment  are  shown  in  Fig.  16,  plotted  as 


the  heat-transfer  parameter,  q^  vRQ/Pg  ,  versus  the  flight  \alocity, 

s  e 

U  .  Here  d  is  the  heating  rate,  R  is  the  nose  radius,  and  P  is  the 
s  e 

stagnation  pressure  at  the  nose  of  the  vehicle.  Also  shown  in  Fig.  16 
are  experimental  data  obtained  in  a  shock  tube,  as  reported  by  Rose  and 
Stark  (Ref.  53). 


The  theoretical  prediction  of  Fay  and  Riddell,  obtained  from  Eq.  40, 
is  shown  in  Fig.  16.  For  comparison  the  correlation  presented  by  Detra, 
Kemp,  and  Riddell  is  also  shown  in  Fig.  16. 

It  is  apparent  from  Fig.  16  that  the  expression  developed  by  Fay  and 
Riddell  can  be  used  with  confidence  down  to  velocities  of  the  order  of 
2,000  fps  if  the  wall  enthalpy  is  properly  considered. 

Table  3  presents  tabulated  values  of  the  experimental  results  obtained 
by  Buglia  (Ref.  50)  for  the  stagnation-point  heat-transfer  rates  for  the 
flight  Mach  numbers  of  2.32,  2.47,  2.63,  2.80,  2.97,  and  3.14.  Also  pre¬ 
sented  are  values  of  the  stagnation-point  heat- transfer  rates  calculated 
from  first,  Fay  and  Riddell  according  to  Eq.  40;  second,  from  Detra,  Kemp, 
and  Riddell,  according  to  Eq.  42;  and  finally,  from  Cohen,  according  to 
Eq.  41. 


For  this  limited  range  of  results,  the  Fay  and  Riddell  equation 
indicates  an  average  error  of  13%  high,  whereas,  Eq.  41  from  Cohen 
indicates  an  average  error  of  10%. 

For  conditions  where  the  enthalpy  at  the  wall  temperature  becomes 
of  the  same  order  as  the  stagnation  enthalpy,  expe.  cental  results  pre¬ 
sented  by  Chauvin  (Ref.  51)  and  Chauvin  and  Maloney  (Ref.  52)  are  of 
interest.  The  heat  transfer  to  a  4-inch  hemisphere  was  measured  for 
several  flow  Mach  numbers  and  stagnation  temperatures.  The  stagnation- 
point  heat-transfer  rates  are  presented  in  Fig.  17  along  with  the  data 
of  Buglia  (Ref.  50)  as  functions  of  the  enthalpy  drop  across  the  boundary 
layer  at  the  stagnation  point.  These  results  may  be  compared  with  the 
theoretical  results  of  Fay  and  Riddell  for  a  Lewis  number  of  unity,  as 
shown  in  Fig.  17. 


50 


NWC  TP  4570 


FVI6HT  VELOCITY,  U.  x  Id'5,  fp» 


FIG.  16.  Stagnation-Point 
Heat  Transfer  for  Air. 


It  thus  becomes  apparent  that  additional  blunt-body  stagnation-point 
heat-transfer  data  are  required  for  flight  Mach  numbers  in  the  range  of 
0.7  to  2.0  in  order  to  verify  the  applicability  of  the  analytical  methods. 


LAMINAR  HEAT  TRANSFER  OVER  SPHERICAL  SECTIONS 

The  distribution  of  heat  transfer  over  the  surface  of  the  hemisphere 
in  the  region  of  the  stagnation  point  and  along  the  surface  of  the  cone 
section  of  the  vehicle  examined  by  Buglia  can  be  compared  in  the  following 
way. 
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Fay  and  Riddell,  and  Cohen  found  that  the  variation  in  fluid 
pr  ertles  between  wall  conditions  and  stagnation  conditions  at  the 
outer  edge  of  the  boundary  layer  must  be  taken  into  account  in  computing 
the  heat-transfer  distribution.  One  method  for  accounting  for  variable 
fluid  properties  has  been  advanced  by  Eckert  (see  Ref.  54  for  a  thorough 
discussion  and  an  extensive  bibliography  of  early  work)  through  the 
introduction  of  the  reference  temperature,  defined 

T*  -  T  +  0.50  /T  -  T  \  +  0.22  /T 

e  (  w  e  /  (aw 

where  the  adiabatic  wall  temperature  is  given  by 

T  -  T  (l  +  :L~-  M  2  r*) 
aw  e  \  2  e  / 

if 

the  recovery  factor,  r  -  Pr  ,  for  the  laminar  boundary  layer,  and 
*  1/3 

r  -  Pr  for  the  turbulent  boundary  layer. 

The  laminar  heat-transfer  distribution  over  the  surface  of  the 
hemisphere,  as  predicted  by  the  various  methods  of  analysis,  are  compared 
on  the  basis  of  the  local  Stanton  number,  defined  as 

°H  '  Vlp‘ue(hav  -  0]  (45) 


as 

-  Te  )  (43) 


(44) 


where  the  reference  density,  p  ,  is  computed  through  the  perfect  gas 
equation  using  the  static  pressure  at  the  outer  edge  of  the  boundary 
layer  and  the  local  reference  temperature. 

The  methods  of  analysis  for  predicting  the  heat-transfer  distribution 
over  the  hemisphere  surface,  which  appear  to  be  the  most  straightforward 
and  yet  within  acceptable  agreement  with  the  experimental  data,  are  the 
methods  of  Korobkin  (Ref.  47)  and  Lees  (Ref.  39). 

From  Sibulkin's  equation  for  the  stagnation-point  heat  transfer, 

Eq.  33,  Korobkin  replaced  the  length,  L,  by  the  wetted  distance  from 
the  stagnation  point,  s,  and  the  velocity  gradient,  6,  by  u  /s  and 
obtained  the  equation 


(46) 
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where  the  properties  are  evaluated  at  the  outer  edge  of  the  boundary 
layer.  This  expression  may  also  be  written  in  terms  of  the  Stanton 
number  by  dividing  through  by  Pr  Res  where  the  Reynolds  number ,  based 
on  the  conditions  at  the  edge  of  the  boundary  layer,  is  written 


Re 

s 


pus 
e  e 


(47) 


Korobkin  thus  obtained  the  expression 


CH 


0.763 


/  -0.6\ 

/  \ 

(Pr  ) 

K) 

(48) 


for  the  Stanton  number  in  the  stagnation  region.  To  account  for  the 
effect  of  wall  temperature  on  the  gas  properties,  this  expression  may 
be  written  in  terms  of  properties  evaluated  at  the  ref'srence  enthalpy: 


4  ■  °-763  (p?)-0,6  K) 


(49) 


for  the  stagnation  region  of  a  hemisphere.  This  relationship  is  conserva¬ 
tive  in  the  computation  of  the  heat-transfer  distribution  over  the  hemis¬ 
phere  surface. 

Lees  (Ref.  39)  presents  the  following  expression  for  the  ratio  of 
the  rate  of  heat  transfer  to  a  given  point  on  the  hemisphere  surface  to 
the  heat  transfer  at  the  stagnation  point: 


% 

% 


26  sin  6  (1  -  1/YM^2)  cos26  +  1/ymJ! 

fO)]*5 


(50) 


54 


-*  r-tgt-  -,-r 


NWC  TP  4570 


3 

A 

t 


where 


f(0) 


6  sin  49  1  -  cos  46 

2  8 


i 

(51) 

I 


+  — —  (e!  -  0  sin  20  +  1  ~  c°s  29  ) 
yM  2  '  2  ' 

Here  0  ie  the  angle  betweer  the  flight  velocity  vector  and  the  hemisphere 
radius  vector.  This  expression  is  applicable  only  in  a  region  of  laminar 
boundary  layer  on  the  hemisphere  surface. 


j 

I 


l 

J 


Thus,  for  the  stagnation  region  of  a  blunt-body  vehicle,  a  comparison 
is  presented  of  the  equation  of  Korobkin,  Eq.  49,  based  on  reference 
properties,  and  the  equation  of  Lees,  Eq.,  50,  with  appropriate  experimental 
results . 


For  blunt-body  vehicles  that  include  either  conical  sections  or 
cylindrical  sections,  the  flat-plate  formula  presented  by  Van  Driest 
(Ref.  55)  should  be  available.  The  flat-plate  heat-transfer  equation 
for  laminar  flow  is  written 


0.332 


(p?) 


i -2/3  /»: 


* 

.  Re 

\  B) 


(52) 


where  the  properties  are  determined  at  the  reference  temperature  and  the 
static  pressure  at  the  outer  edge  of  the  boundary  layer. 

The  experimental  data  obtained  by  Buglla  (Ref.  50)  for  the  heat- 
transfer  distribution  over  the  hemisphere  and  cone  surfaces  were  reported 
in  the  form 


s* 

CH  "  pu/h  -  h  \ 

00  00  col  a  w  I 


and  were  corrected  to  the  reference  temperature  base  by 
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The  experimental  values  obtained  by  Buglia  (Ref.  50)  are  presented  in 
Fig.1 13(a)  through  18(f)  as  a  function  of  the  Reynolds  number,  defined  as 


,  * 

p  u  u 

e 


(55) 


where  s  is  the  wettecl  length  from  the  stagnation  point:  to  the  point  in 
question.  The  static  pressure  at  the  edge  of  the  boundary  layer,  employed 
in  the  calculatior  of  the  local  temperature  and  velocity,  wcs  obtained 
from  the  experimental  pressure  distribut  ion.  Isentropic  expansion  from 
the  stagnation  point  to  the  local  pressure  was  then  assumed  tc  calculate 
the  local  temperature  and  velocity.  The  data  points  presented  in  Fig.  18(a) 
through  18(f)  on  the  cone  section  are  indicated  with  a  flag.  It  is  inter¬ 
esting  to  note  in  the  vicinity  of  the  juncture  between  hemisphere  and 
conical  sections  that  the  Reynolds  number,  based  on  reference  conditions, 
decreases.  As  the  flow  is  compressed  in  the  region  of  the  conical  section, 
the  boundary  layer  undergoes  transition  from  laminar  to  turbulent  flow, 
thus  showing  a  significant  increase  in  the  heat-transfer  rate. 

The  theoretical  predictions  of  Korobkin,  Eq.  49:  Lees,  Eq.  50;  and 
Van  Driest,  Eq.  52,  are  also  included  in  each  of  Fig.  18(a)  through  18(f) 
for  comparison  with  the  experimental  results. 

The  results  presented  in  Fig.  18(a)  through  18(f)  indicate  that  the 
Lee6  relationship,  based  on  reference  values,  more  closely  predicts  the 
experimental  data  for  laminar  flow  over  both  the  hemisphere  and  the  initial 
portion  of  the  cone  section  of  the  body.  Also,  for  a  case  where  laminar 
flow  persists  over  a  significant  portion  of  the  conical  sectio.i,  it  is 
apparent  that  the  Van  Driest  equation,  Eq .  52,  may  be  used  to  predict  the 
heat  transfer. 


It  is  also  apparent  that  fox  the  Mach  number  range  from  2  to  3,  if 
laminar  flow  persists  over  the  entire  hemisphere  section,  Lees'  expression 
will  account  adequately  for  the  pressure  gradient  effect  over  the  juncture 
between  the  hemisphere  and  the  conical  portion  of  the  body. 

The  numerical  computations  corresponding  to  the  results  prese  ited  in 
Fig.  18(a)  through  18(f)  are  given  in  Tables  4(a)  through  4(f).  Each  table 
contains  the  numerical  results  presented  in  the  corresponding  figure.  The 
parameters  of  interest,  for  example,  the  local  Reynolds  number  based  on 
wetted  length  and  reference  conditions,  are  presented  in  terms  of  the  non- 
dimensional  wetted  distance  from  the  stagnation  point  for  each  Mach  number. 
The  results  presented  in  these  tables  indicate  that  the  best  agreement 
with  the  experimental  Stanton  number  distribution  is  achieved  by  using  the 
Fay  and  Riddell  stagnation-point  heat-transfer  equation,  Eq.  40,  with  Lees' 
expression  for  the  distribution  of  the  local  heat-transfer  equation,  Eq.  50, 
and  the  density  of  the  air  computed  with  the  value  of  the  reference  tempera¬ 
ture  and  local  static  pressure.  This  method  will  be  presented  in  detail 
in  a  later  section. 
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FIG.  18.  Heat-Transfer  Distribution  Over  a  Smooth 
Hemispherically  Capped  Cone  for  Flight  Mach  Numbers 
Buglia,  Ref.  50. 
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TABLE  4(a).  Experimental  and  Theoreti:al  Heat-Transfer 
Distributions  Over  a  Hemispherically  Capped  Cone 

M  -  2.32,  P  -  11.0  psia,  T  -  500. 8°R 

00  7  oo  r  1  00 


s/2  Ro 

R* 

es 

x  10'6 

C*  (1) 

GXp 

x  103 

CH  (2) 

hl 

x  103 

cHl 

Le 

x  103 

V4’ 

x  103 

cHi  <5’ 
Lp-E 
x  103 

0.090 

0.479 

1.15 

1.242 

1.458 

1.155 

1.355 

0.150 

1  195 

0.67 

0.805 

0  942 

0.748 

0.876 

0.210 

2.162 

0.37 

0.592 

0.686 

0.550 

0.636 

0.270 

3.411 

0.40 

0.440 

0.502 

C .  40° 

0.466 

0.330 

4 .  482 

0.32 

0.37? 

0.419 

0.351 

0.389 

0.390 

5.323 

0.32 

0.335 

0.358 

0.312 

0.333 

0.520 

5.696 

0.20 

0.275 

0.267 

0.256 

0.249 

0. 660 

3.056 

0.28 

0.329 

0.264 

0.306 

0.246 

0.810 

3.304 

0.50 

0.332 

0.259 

0.309 

0.241 

1.050 

5.758 

0.22 

0.186 

0.153 

0.173 

0.142 

1.400 

7.415 

1.57 

0.215 

0.170 

0.200 

0.158 

Experimental  Stanton  Number,  Cj.  ■  — jj - - -  ;  Buglia, 

exp  p  u  |h  -  h  ) 
e \  aw  w / 

Ref.  50. 

(2' 

'Stanton  Number  based  on  the  Fay  and  Riddell  stagnation-point 
heat  transfer,  Lees'  distribution  of  local  heat  transfer,  and  reference 
conditions . 

(3) 

Stanton  Number  based  on  the  Fay  and  Riddell  stagnation-point  heat 
transfer,  Lees'  distribution  of  local  heat  transfer,  and  conditions  at 
the  edge  of  the  boundary  layer. 

(4) 

Stanton  Number  based  on  the  experimental  stagnation-puint  heat 
transfer,  Lees'  distribution  of  local  heat  transfer,  and  reference 
conditions . 

^ Stanton  Number  based  on  the  experimental  stagnation-point  heat 
transfer,  Lees'  distribution  of  local  heat  transfer,  and  conditions 
at  the  outer  edge  of  the  boundary  layer. 


*  i 
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TABLE  4(b).  Experimental  end  Theoretical  Heat -Transfer 
Distributions  Over  a  Hemispherically  Capped  Cone 

M  -  2.47,  P  -  1C. 8  psia,  T  «  500°R 
00  '00  *  '  00 


s/2  R 

o 

R* 

es 

x  10-6 

CH  (1) 
exp 

x  103 

C£  (2) 
kl 

x  103 

C„L  (3) 
x  lo3 

(4) 

Lre 

x  103 

Cfi.  (5) 

l*e 

x  103 

0.090 

0.499 

1.09 

1.221 

1.450 

1.132 

1.345 

0.150 

1.245 

0.62 

0.794 

0.939 

0.736 

0.871 

0.210 

2.252 

0.40 

0.585 

0.685 

0.542 

0.635 

0.270 

3.552 

0.41 

0.437 

0.505 

0.406 

0.468 

0.330 

4.800 

0.25 

0.344 

0.393 

0.319 

0.364 

0.390 

5.852 

0.29 

0.285 

0.318 

0.264 

0.295 

0.520 

6.216 

0.18 

0.250 

0.252 

0.232 

0.234 

0.660 

3.129 

0.47 

0.368 

0.297 

0.342 

0.275 

0.810 

3.404 

0.81 

0.375 

0.296 

0.348 

0.274 

1  050 

5.944 

0.33 

0.191 

0.158 

0.177 

0.147 

.1.400 

7.655 

1.55 

0.220 

0.176 

0.204 

0.163 

Cl)  ^  6XD 

Experimental  Stanton  Number,  *  — - *- -  ;  Buglia, 

exp  p  u  (h  -  h  ) 
e  \  aw  w  / 

Ref.  50. 

(2) 

Stanton  Number  based  on  the  Fay  and  Riddell  stagnation-point 
heat  transfer,  Lees'  distribution  of  local  heat  transfer,  and  reference 
conditions . 

(3) 

Stanton  Number  based  on  the  Fay  and  Riddell  stagnation-point 
heat  transfer,  Lees'  distribution  of  local  heat  transfer,  and  conditions 
at  the  edge  of  the  boundary  layer. 

(4) 

Stanton  Number  based  on  the  experimental  stagnation-point  beat 
transfer,  Lees'  distribution  of  local  heat  transfer,  and  reference 
conditions . 

^Stanton  Number  based  on  the  experimental  stagnation-point  heat 
transfer,  Lees'  distribution  of  local  heat  transfer,  and  conditions 
at  the  outer  edge  of  the  boundary  layer. 
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TABLE  4(c) .  Experimental  and  Theoretical  Heat-Transfer 
Distributions  Over  a  Hemispherically  Capped  Cone 

M  -  2.63,  P  -  10.6  psia,  T  -  499°R 

00  OO  *  00 


s/2  R 

0 

R* 

ct 

x  10-A 

ca  (l) 

exp 
x  103 

CS  (2) 

x  103 

CHL  O) 

Le 

x  It- 

A  (4> 

LRE 
x  103 

CHL  (5) 
LeE 
x  103 

0.090 

0.515 

0.97 

1.204 

1.439 

1.076 

1.285 

0.150 

1.271 

0.61 

0.810 

0.959 

0.724 

0.857 

0.210 

2.300 

0.39 

0.596 

0.698 

0.532 

0.623 

0.270 

3.626 

0.40 

0.445 

0.513 

0.397 

0.458 

0.330 

4.898 

0.27 

0.348 

0.397 

0.311 

0.355 

0.390 

5.969 

0.30 

0.286 

0.320 

0.256 

0.286 

0.520 

6.334 

0.32 

0.243 

0.246 

0.217 

0.219 

0.660 

3.342 

1.05 

C.302 

0.250 

0.270 

0.224 

0.810 

3.554 

1.30 

0.308 

0.246 

0.275 

0.220 

1.050 

6.088 

0.96 

0.193 

0.160 

0.172 

0.143 

1.400 

7.844 

1.53 

0.222 

0.178 

0.198 

0.159 

(1)  ^  BXP 

Experimental  Stanton  Number,  C^  »  — j c - -  ;  Buglia, 

exp  p  u  /h  -  h 
v  v  e\  aw  w/ 

Ref.  50. 

(2) 

Stanton  Number  based  on  the  Fay  and  Riddell  stagnation-point 
heat  transfer,  Lees'  distribution  of  local  heat  transfer,  and  reference 
conditions . 

(3) 

Stanton  Number  based  on  the  Fay  and  Riddell  stagnation-point 
heat  transfer,  Lees'  distribution  of  local  heat  transfer,  and  conditions 
at  the  edge  of  the  boundary  layer. 

(A) 

Stanton  Number  based  on  the  experimental  stagnation-point  heat 
transfer,  Lees'  distribution  of  local  heat  transfer  and  reference 
conditions . 

^Stanton  Number  based  cn  the  experimental  stagnation-point  heat 
transfer,  Lees'  distribution  of  local  heat  transfer,  and  conditions  at 
the  outer  edge  of  the  boundary  layer. 
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TABLE  4(d).  Experimental  and  Theoretical  Heat-Transfer 
Distributions  Over  a  Hemispherically  Capped  Cone 

M  »  2.80,  P  -  10.5  psia,  T  -  498°R 

00  '00  ‘  00 


s/2  R 

0 

R* 

es 

x  10“6 

c*  (1) 
exp 
x  103 

CH  (2) 

HL 

x  103 

CH,  (3) 

Le 

x  103 

ejS.  (4) 

^lre 

x  103 

CHL  <5> 
LeE 

x  103 

0.090 

0.563 

0.82 

1.150 

1.434 

0.997 

1.244 

0.150 

1.404 

0.46 

0.752 

0.934 

0.652 

0.810 

0.210 

2.565 

0.43 

0.540 

0.667 

0.468 

0.578 

0.270 

4.038 

0.36 

0.408 

0.496 

0.354 

0.430 

0.330 

5.287 

0.22 

0.352 

0.416 

0.305 

0.360 

0.390 

6.191 

0.42 

0.324 

0.366 

0.281 

0.317 

0.520 

6.496 

0.48 

0.280 

0.284 

0.243 

0.246 

0.660 

3.412 

1.94 

0.336 

0.280 

0.292 

0.243 

0.810 

3.703 

1.95 

0.325 

0.264 

0.281 

0.229 

1.050 

6.500 

1.52 

0.210 

0.180 

0.182 

0.156 

1.400 

8.362 

1.50 

0.239 

0.197 

0.207 

0.171 

( 1)  *  pyp 

Experimental  Stanton  Number,  *»  — j K -  :  Buglia, 

exp  p  u  /h  -  h  ) 
e\  aw  w f 

Ref.  50. 

(2) 

Stanton  Number  based  on  the  Fay  and  Riddell  stagnation-point 
heat  transfer,  Lees'  distribution  of  local  heat  transfer,  and  reference 
conditions . 


(3) 

Stanton  Number  based  on  the  Fay  and  Riddell  stagnation-point 
heat  transfer,  Lees'  distribution  of  local  heat  transfer,  and  conditions 
at  the  edge  of  the  boundary  layer. 

(4) 

Stanton  Number  based  on  the  experimental  stagnation-point  heat 
transfer,  Lees'  distribution  of  local  hear  transfer,  and  reference 
conditions . 

Stanton  Number  based  on  the  experimental  stagnation-point  heat 
transfer,  Lees'  distribution  of  local  heat  transfer,  and  conditions  at 
the  outer  edge  of  the  boundary  layer. 
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TABLE  4(e) .  Experimental  and  Theoretical  Heat-Transfer 
Distributions  Over  a  Hemispherically  Capped  Cone 

M  -  2.97,  Pm  -  10.2  psia,  T  -  497°R 
00  '00  ‘  00 


b/2  R 

0 

R* 

es 

*  io“6 

CH  (1) 

eXP  i 

x  103 

«  (2) 
h 

x  103 

C«L  (3) 

Le 

x  103 

V4’ 

X  103 

CV5’ 

x  103 

0.090 

0.581 

0.74 

1.134 

1 

1.432 

0.972 

1.227 

0.150 

1.433 

0.45 

0.764 

0  955 

0.655 

0.818 

0.210 

2.589 

0.43 

0.564 

0.698 

0.484 

0.598 

0.270 

4.117 

0.36 

0.414 

0.507 

0.355 

0.435 

0.330 

5.497 

0.22 

0.339 

0.408 

0.291 

0.350 

0.390 

6.552 

0.71 

0.298 

0.346 

0.255 

0.296 

0.520 

6.958 

0.94 

0.249 

0.262 

0.214 

0.225 

0.660 

3.428 

2.36 

0.332 

0.277 

0.285 

0.237 

0.810 

3.720 

2.36 

0.314 

0.256 

0.270 

0.219 

1.050 

6.585 

1.69 

0.217 

0,187 

0.186 

0.160 

1.400 

8.535 

1.62 

0.243 

0.202 

0.208 

0.173 

( 1)  ^  gxp 

Experimental  Stanton  Number,  C^  -  — - c -  ;  Buglia, 

exp  p  u  (h  -  h  ) 
e\  aw  w/ 

Ref.  50. 

(2) 

Stanton  Number  based  on  the  Fay  and  Riddell  stagnation-point 
heat  transfer,  Lees  distribution  of  local  heat  transfer,  and  reference 
conditions . 

(3) 

Stanton  Number  based  on  the  Fay  and  Riddell  stagnation-point 
heat  transfer,  Lees'  distribution  of  local  heat  transfer,  and  conditions 
at  the  edge  of  the  boundary  layer. 

(4) 

Stanton  Number  based  on  the  experimental  stagnation-point  heat 
transfer,  Lees'  distribution  of  local  heat  transfer,  and  reference 
conditions . 

^ Stanton  Number  based  on  the  experimental  stagnation-point  heat 
transfer,  Lees'  distribution  of  local  heat  transfer,  and  conditions  at 
the  outer  edge  of  the  boundary  layer. 
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TABLE  4(f).  Experimental  and  Theoretical  Heat-Transfer 
Distributions  Over  a  Hemispherically  Capped  Cone 

M  -  3.14,  -  10.0  psia,  T  -  495°R 

00  '00  *  00 


c/2  R 

o 

R* 

es 

*  10-6 

Cg  (1) 
exp 
x  103 

CH  (2) 
x  103 

cHi  (3) 

Le 

x  103 

V4’ 

x  103 

cHi  (5) 

l«e 

x  103 

0.090 

0.598 

0.63 

1.118 

1.429 

0.925 

1.183 

0.150 

1.491 

0.45 

0.732 

0.933 

0.606 

0.772 

0.210 

2.693 

0.42 

0.542 

0.632 

0.448 

0.565 

0.270 

4.194 

0.35 

0.421 

0.518 

0.348 

0.429 

0.330 

5.597 

0.21 

0.344 

0.417 

0.285 

0.345 

0.390 

6.867 

1.04 

0.281 

0.335 

0.233 

0.277 

0.520 

6.766 

1.344 

0,278 

0.287 

0.230 

0.237 

0.660 

3.419 

2.29 

0.334 

0.277 

0.276 

0.229 

0.810 

3.683 

2.52 

0.313 

0.252 

0.259 

0.209 

1.050 

6.476 

1.66 

0  ’36 

0.201 

0.196 

0.166 

1.400 

8.584 

1.40 

0.254 

0.211 

0.210 

0.175 

(1)  *  CXP 

Experimental  Stanton  Number,  CR  ■  — j - c - -  ;  Buglla, 

exp  p  u  fh  -  h  ) 
e(  aw  w/ 

Ref.  50. 

(2) 

Stanton  Number  based  on  the  Fay  and  Riddell  stagnation-point 
heat  transfer,  Lees'  distribution  of  local  heat  transfer,  and  reference 
conditions . 

(3) 

Stanton  Number  based  on  the  Fay  and  Riddell  stagnation-point 
heat  transfer,  Lees'  distribution  of  local  heat  transfer,  and  conditions 
at  the  edge  of  the  boundary  layer. 

(4) 

Stanton  Number  based  on  the  experimental  stagnation-point  heat 
transfer,  Lees'  distribution  of  local  heat  transfer,  and  reference 
conditions . 

^Stanton  Number  based  on  the  experimental  stagnation-point  heat 
transfer,  Loes'  distribution  of  local  heat  transfer,  and  conditions  at 
the  outer  edge  of  the  boundary  layer. 
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It  should  be  pointed  out  that  the  experimental  results  used  for  the 
comparisons  made  above  were  obtained  on  a  smooth  (roughness  from  2  to 
5  y-inches)  hemispherically  capped  cone.  The  polished  surface  condition 
apparently  delayed  transition  from  laminar  to  turbulent  flow  for  this 
particular  vehicle;  hence,  the  Reynolds  numbers  for  which  laminar  flow 
persists  are  high,  on  the  order  of  5  x  10^.  These  data  are  thus  presented 
only  to  compare  the  various  methods  of  predicting  laminar  blunt-body 
heat-transfer  distributions  with  available  experimental  results. 

Additional  experimental  data  for  the  heat-transfer  distribution  over 
hemispherically  capped  cylinders  have  been  reported  by  Chauvin  and  Maloney 
(Ref.  52)  and  Korobkin  (Ref.  47).  However,  an  insufficient  amount  of 
supporting  test  information  prohibited  comparison  of  these  data  with  the 
methods  previously  discussed.  These  data  are,  therefore,  presented  in 
Fig.  19  as  the  experimental  Stanton  number  based  on  local  conditions  at 
the  edge  of  the  boundary  layer  versus  the  local  Reynolds  number.  The 
Mach  number  range  covered  by  the  wind  tunnel  tests  of  Chauvin  and  Maloney 
was  from  1.62  to  3.04,  rhile  Korobkin  obtained  a  Mach  number  range  from 
1.90  to  4.87.  These  data  Indicate  a  fairly  uniform  distribution  about 
a  line  predicted  by  Korobkin 's  equation,  Eq.  49.  It  should  be  noted 
that  in  the  tests  reported  by  Chauvin  and  Maloney,  and  Korobkin,  the 
boundary  layer  underwent  transition  from  laminar  to  turbulent  at  a  local 
Reynolds  number  of  approximately  1  x  10^,  based  upon  the  wetted  length 
from  the  stagnation  point  and  the  properties  evaluated  at  local  conditions 
at  the  edge  of  the  boundary- layer  heat-transfer . 


FIG.  19.  Local  Stanton  Number  Distribution  Over  a  Hemispherically 
Capped  Cylinder  Obtained  in  Wind-Tunnel  Tests;  Chauvin  and  Maloney, 
Ref.  52,  and  Korobkin,  Ref.  47. 
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Cresci,  MacKenzie,  and  Libby  (Ref.  57)  have  compared  laminar  heat- 
transfer  data  with  the  theory  of  Lees  (Ref.  39)  by  using  the  stagnation- 
point  value  of  Fay,  Riddell,  and  Kemp  (Ref.  56)  for  Mach  numbers  above  3. 
The  data  are  presented  in  Ref.  57  in  terms  of  a  Nusselt  number  defined  by 


Nu  “ 


(56) 


and  a  Reynolds  number  parameter 


Re  =  Re 


v/p  /p  h 

s  s  s 
e  e 


Re  P' 
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In  terms  of  these  parameters,  the  heat-transfer  distribution  is  represented 
by 
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Corresponding  to  this  equation,  the  Reynolds  number  based  on  momentum 
thickness  is  given  by  the  theory  of  Lees : 


where 


and 


Re, 


p  u  0 
e  e 


e  -  f“  (PM/Peliej  [l  -  (u/ue)  ]  Ay 


(61) 


(62) 


(63) 


Comparison  of  Eq.  55  and  experimental  data  is  presented  in  Ref.  57 
with  excellent  agreement.  Since  these  experimental  results  corresponded 
to  a  flight  Mach  number  of  approximately  20,  they  are  not  presented  here. 


HEAT  TRANSFER  IN  THE  TRANSITION  REGION 


An  extensive  discussion  of  boundary-layer  transition  is  presented 
by  Wilson  in  Ref.  58.  This  work,  however,  is  largely  concerned  with 
transition  phenomenon  on  pointed  bodies  of  revolution.  In  fact,  the 
observation  is  made  in  this  reference  that  few  data  are  available  for 
blunt-body  heat  transfer  in  the  transition  region.  Despite  this,  a 
computational  procedure  has  been  compared  with  experimental  results  for 
the  heat  transfer  in  the  transition  regime  by  Cresci,  Mackenzie,  and 
Libby  (Ref.  57),  by  Nardo  and  Sadler  (Ref.  59),  and  by  Stetson  (Ref.  60). 
In  anticipation  of  its  use  in  transition  flow,  the  momentum  equation  for 
the  axisymmetric  boundary  layer  can  be  written  in  the  form 


dRe, 


j/ds  -  ^Re/agl5^pu/ij^Cf/2j  -  Regd(ln  r  y)/ds 


(64) 


68 


NWC  TP  4570 


where  C,  -  2t  /p  u  2,  the  local  skin-friction  coefficient  a  •  P  /p  h 
fw  ee  sees 

e  e  e 

The  analysis  of  Ref.  57  utilizes,  for  fully  developed  turbulent  flows, 
the  extension  of  the  local  Prandtl  incompressible  skin-friction  law  for 
the  incompressible  "cold-wall"  case  in  the  form 


Cf/2  -  y(0.013)  (65) 


Numerical  integration  of  Eq.  64  along  with  a  skin-friction  law  equivalent 
to  Eq.  65  is  required.  In  accordance  with  the  suggestion  of  Persch 
(Ref.  61)  and  including  the  modification  suggested  by  Ref.  55,  the  skin- 
friction  becomes 


Cf/2 


_L- 

U 10.013  Refi  -  JI  Re 


-2 


0 


(66) 


The  constant  II  in  this  expression  is  equivalent  from  the  laminar 
flow  such  that  the  friction  coefficient,  Cf/2,  remains  continuous  at 
the  transition  point  (Ref.  57), 


n 


(67) 


If  Eq.  64  and  66  are  combined,  and  Eq.  64  is  integrated  from  the 
transition  point,  then  the  skin-friction  coefficient  and  the  momentum- 
thickness  distribution  are  known.  When  Reynolds  analogy  is  modified 
to  account  for  nonunity  Prandtl  number  effects,  then  the  heat-transfer 
distribution  may  be  found  from 


4* 


Pr 


-(2/3) 


(68) 


Comparisons  of  this  method  of  analysis  with  experimental  wind- 
tunnel  data  has  been  made  by  Cresci,  Mackenzie ,  and  Libby  (Ref.  57) 
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lind  by  Nardo  and  Sadler  (Rtf.  59).  These  comparisons  were  based  upon 
3:he  Nusselt  number  evaluated  from  the  expression 


Figures  20(a)  through  20(e)  show  the  variation  of  Nusselt  number  along 
a  meridian  of  the  hemisphere  used  in  the  wind-tunnel  tests .  These 
results  were  obtained  for  Re  from  0.97  x  10^  to  5.1  x  10^  with  the 
total  temperature  variation  from  approximately  1200— 1900°R.  The  tunnel 
Mach  numbers  were  not  given.  These  results  indicate  a  qualitative 
agreement  of  the  computational  procedure  with  the  experimental  results. 
It  is  readily  apparent,  however,  that  appropriate  free-f light  transition 
heat-transfer  data  will  be  required  before  the  method  can  be  used  with 
confidence  for  design  purposes.  No  such  flight  data  are  presently 
available  in  the  open  literature. 


0  0.2  O.h  0.6  0.8  1.0  1.2  I.*  1.6 
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FIG.  20(a) .  Variation  of  Nusselt  Number  Along 
Meridian  of  Hemisphere,  Re  **  0.97  x  10^;  Nardo 
and  Sadler,  Ref.  59. 
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FIG.  20.  Variation  of  Nusselt  Number  Along 
Meridian  of  a  Hemisphere;  Nardo  and  Sadler, 
Ref.  59. 


TURBULENT  HEAT  TRANSFER  ON  BLUNT  BODIES  OF  REVOLUTION 
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Fully  developed,  turbulent  he&t  transfer  over  blunt -body  vehicles 
has  been  compared  to  the  prediction  of  the  flat -plate  reference-enthalpv 
method  by  Cresci,  MacKenzie,  and  Libby  (Ref.  57)  for  wind-tunnel  results 
and  by  Zoby  and  Sullivan  (Ref.  62)  ,  and  by  Eckert  (Ref.  63)  for  the 
results  of  free-f light  tests. 

The  correlations ,  presented  by  these  authors ,  vrcre  made  by  relating 
the  local  skin-friction  coefficient  to  the  heating  rate  by  a  modified 
Reynolds  analogy  expressed  as 


f 


(70) 


To  evaluate  the  Stanton  number,  C^,  for  turbulent  heat  transfer,  the 
local  skin-friction  coefficient,  Cf,  must  be  pioperly  calculated.  The 
Blasius  relation  fcr  the  incompressible  skin-friction  coefficient  is 


-  0.0296  (Re  )“0-2  (71) 


for  the  Reynolds  number  range  105  <_  Res  £  109,  and  the  Schultz-Grunow 
incompressible  skin-friction  relation  is 


(-§),  •  °'185  (lo*io  K's)'2-584 


for  the  Reynolds  number  range  105  £  Res  £  109.  Zoby  and  Sullivan  (Ref  62) 
point  out  that  both  relationships  correlate  the  data  below  a  Reynolds 
number,  Res,  of  1G?  but  that  the  Schultz-Grunow  relation  applies  over  the 
entire  range  of  Reynolds  numbers  from  105  to  109. 

Compressibility  effects  were  accounted  for  by  again  evaluating  the 
flow  properties  at  reference  conditions.  The  local  skin-friction  coeffi¬ 
cient  is  then  written 


cf/2  -  0.0296  (p*/Pe)°*8  (y*/yej0,2 
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and 


-2.584 


(74) 


V2  -  0.185  (p*/pe)  (log10  R 


for  the  Blaslus  and  Schultz-Grunow  relations,  respectively. 

By  substituting  these  expressions  into  Eq.  70,  the  following 
expressions  for  the  heating  rate  are  obtained: 


V  -  0.0296  peuc(hm  -  hu)(Pr*)-2/3(Res)-°-2(p*/pe)0'8(p*/pe)0-2  (75) 


and 


S»  -  0.185  U,u,(hw  -  ’>w)(Pr*)'2/3(lo810  »'*)  2‘584(p‘/pJ  <») 


It  is  easily  shown  that  Eq.  70  may  be  written  in  terms  of  the  Eckert 
flat-plate  reference-enthalpy  formula.  If 
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Re 
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pus 
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(77) 


then 
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CH  "  * 


u  /  h  -  h  \ 
e \  sw  w/ 


-  0.0296  (pr)~2/3  ^ Reg^-0’ 2 


(78) 


The  Reynolds  number  is  based  upon  properties  evaluated  at  reference 
conditions,  the  velocity  at  the  outer  edge  of  the  boundary  layer,  and 
the  corresponding  flat-plate  wetted  length. 
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Experimental  results  for  turbulent  heat  transfer  to  blunt  bodies 
of  revolution  obtained  In  wind  tunnel  tests,  presented  In  Ref.  57,  are 
shown  In  Fig.  21.  These  data  were  obtained  from  the  shrouded-model 
wind-tunnel  technique  for  tunnel  total  temperature  of  1200— 1900°R, 
hence,  tunnel  Mach  numbers  were  not  reported. 

Shown  also  in  Fig.  21  are  results  obtained  by  Buglla  (Ref.  50)  from 
free-f light  tests  of  the  smooth  hemispherlcally  capped  cone  over  a  Mach 
number  range  from  2.32  to  3.14. 


FIG.  21.  Variation  of  Local  Stanton  Number  Based  on 
the  Flat-Plate  Reference-Enthalpy  Method  With  the 
Local  Reynolds  Number;  Cresci,  Mackenzie,  and  Libby, 
Ref.  57;  and  Buglla,  Ref.  50. 


The  correlation  study  of  Zoby  and  Sullivan  (Ref.  62)  Includes  some 
turbulent  heat-transfer  data.  These  data  cover  a  Mach  number  range  from 
2.98  to  5.3  and  include  some  13  data  points.  It  should  be  pointed  out 
that  blunt-body,  turbulent  heat-transfer  data  are  relatively  scarce  as 
compared  to  sharp  cone  results,  and  that  additional  turbulent  heat- 
transfer  results  are  needed  below  a  Mach  number  of  3  to  obtain  more 
confidence  in  the  engineering  heat-transfer  computational  methods. 
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Eckert  (Ref.  63)  has  presented  a  correlation  of  experimental  results 
for  turbulent  flow  over  flat  plates  as  a  function  of  Mach  numbers  in 
the  range  from  0  to  approximately  5.  The  experimental  results  are  not 
pertinent  to  this  report  because  they  were  obtained  for  a  flat  plate, 
but  the  trend  of  the  limited  data  is  important.  Eckert  found  that  the 
ratio  of  the  actual  to  the  constant-property  turbulent-friction  factor 
for  a  plate  at  the  recovery  temperature  scattered  rather  widely  from  the 
theoretical  values  as  the  flow  Mach  number  decreased  from  approximately 
2  to  approximately  unity.  If  this  trend  is  typical  of  heat-transfer 
results  in  this  Mach  number  range,  then  the  method  of  analysis  for  blunt- 
body  heat  transfer  in  the  turbulent  regime  may  be  in  question.  Hence, 
even  though  satisfactory  agreement  for  fully  developed  turbulent  heat- 
transfer  results  may  be  obtained  from  the  flat-plate  reference-enthalp, 
method  for  Mach  numbeic.  greater  than  1.6,  one  must  be  cautious  in  arbi¬ 
trarily  extrapolating  the  results  to  Mach  numbers  close  to  unity.  Clearly, 
additional  experimental  results  are  required  in  this  range. 

For  rapid  calculation  of  turbulent  skin-friction  and  heat- t  ransfer 
results  for  high  Mach  number  vehicles,  Knuth  (Ref.  64),  and  Neal  and 
Bertram  (Ref.  65)  have  presented  turbulent  skin-friction  and  heat-transfer 
charts  adapted  from  the  Spalding  and  Chi  (Ref.  66)  method  for  determining 
the  drag  of  a  compressible,  turbulent  boundary  layer  on  a  smooth  flat 
plate.  This  method  is  presented  for  a  Mach  number  range  from  0  to  20, 
considering  gross  effects;  hence,  its  usefulness  in  a  limited  range  of 
Mach  numbers  from  0.7  to  5  where  detail  is  important  becomes  questionable. 


METHOD  OF  COMPUTATION  OF  THE  HEAT-TRANSFER  DISTRIBUTION 

The  calculation  of  the  stagnation-point  heat-transfer  rate  and  the 
distribution  of  heat-transfer  rate  over  the  surface  of  a  blunt-body 
vehicle  are  preferably  performed  by  a  digital  computer.  The  slide  rule 
can  be  used  for  initial  estimates,  but  the  computations  can  become 
tedious.  The  following  discussion  describes  the  method  of  calculation 
and  the  assumptions  made. 

The  geometrical  configuration  is  assumed  to  be  a  spherically  capped 
frustrum  of  a  cone. 

Effects  of  wall  temperature,  Mach  numbers,  and  variable  fluid 
properties  are  accounted  for  by  evaluating  the  fluid  properties  at  the 
appropriate  reference  conditions.  Correspondingly,  the  theoretical 
variation  of  Stanton  number  versus  Reynolds  number,  based  on  the  reference- 
enthalpy  method,  is  a  straight  line  on  a  log-log  graph,  either  for  laminar 
flow,  Eq.  49,  ov  turbulent  flow,  Eq.  78.  The  state  of  the  boundary  layer 
may  be  readily  ascertained  from  such  a  plot,  as  indicated  in  Fig.  18  (a — f) 
through  20(a— e). 
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To  summarize  the  method,  the  flow  at  the  outer  edge  of  the  boundary 
layer  Is  assumed  to  have  passed  through  a  normal  shock  wave,  followed 
by  lsentrop.ic  expansion  to  the  theoretical  value  of  the  surface  pressure 
found  from  Eq.  12.  Fluid  properties  are  evaluated  at  the  reference 
condition,  thereby  eliminating  the  dependence  of  Stanton  number  and 
Reynolds  number  on  Mach  number,  wall  temperature,  and  variable  fluid 
properties.  Air  properties  are  taken  from  Keenan  and  Kayes  "Gas  Tables," 
whereas  the  viscosity  Is  evaluated  by  using  Sutherland's  equation. 

The  Input  Information  required  for  the  computation  of  heat-transfer 
distributions  is  assumed  to  be  (1)  vehicle  nose  radius,  (2)  flight 
velocity,  (3)  ambient  temperature,  (4)  ambient  density,  and  (5)  wall 
temperature. 

Stagnation-Point  Heat-Transfer  Rate 

The  calculation  of  the  stagnation-point  heat-transfet  -ate  proceeds 
as  follows : 

1.  Compute  the  flight  Mach  number,  1^,  from  the  flight  velocity, 
V^fps)  ,  and  the  ambient  air  temperature,  Tw(#R) : 


2.  Compute  the  pressure  ratio,  P  /P  ,  across  the  normal  shock 

8  00 
e 

wave  from  the  flight  Mach  number  M^: 


|[(y+l)/2]  Ma 


y/y-i 


[2y/(rfl)  Hj  -  (Y-l)/(Y+l)j 


1/Y-l 


For  y  •  1.40,  this  equation  reduces  to 


P 


00 
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3.  Compute  the  ratio  of  the  total  temperature  at  the  outer  edge 
of  the  boundary  layer  at  the  stagnation  point  to  the  free-stream  static 
temperature: 


1  + 


M 


4.  From  the  gas  tables,  extract  the  enthalpy,  haw,  corresponding 
to  the  temperature,  T  (°R),  and  the  enthalpy,  h*,,  corresponding  to 

e 

the  wall  temperature,  Tw  (°R).  For  low  temperatures,  corresponding  to 
<  5,  the  specific  heat  at  constant  pressure  may  be  assumed  constant 
at  a  value  of  Cp  ■  0.246  Btu/lbm°R. 

5.  Compute  the  density  at  the  edge  of  the  boundary  layer  at  the 
stagnation  point,  “  P  /R  Tfl  lbm/ft.3,  and  the  density  of  air  at 

s  e  e 

the  wall  temperature  at  the  nose  of  the  vehicle 


P 

“e  lbm 

Pw  3 

8  RT  ftJ 
ws 


where  R  Is  the  gas  constant  for  air. 

6.  Compute  the  viscosity  of  the  air  at  the  stagnation-point 
temperature,  T  (°R) ,  and  at  the  wall  temperature  at  the  nose  of 
e 

the  vehicle,  T  (°R) ,  from  Sutherland's  formula 
w 

s 


y 


2.270  x  10~8  (T)3/2 
T  198.6 


lb^-sec 

ft2 


where  T  is  the  temperature  in  degrees  Rankine. 
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7.  Compute  the  stagnation-point  heat-transfer  parameter,  ^  v  Ro/Pfl 

a  8e 

from  the  results  of  Fay  and  Riddell,  Eq.  40: 


For  a  given  nose  radius,  R0,  this  result  then  yields  the  stagnation- 
point  heat-transfer  rate. 

Heat-Transfer  Distribution  Over  Hemisphere  Surface 

To  compute  the  heat-transfer  distribution  over  the  hemisphere 
surface,  the  expression  of  Lees,  Eq.  50  and  51  may  be  used  as  follows: 

1.  Compute  the  heat  transfer  on  the  hemispherical  cap  at  the  angle 

0  from  Lees  relationship,  Eq.  50  and  51.  The  angle  0  is  the  angle  between 
the  flight  direction  and  the  radiu3  vector  to  the  surface  station  of 
interest.  The  stagnation-point  heat-transfer  rate  is  obtained  from  the 
previous  calculation. 

An  alternate  method  for  computing  the  local  value  of  the  heat  trans¬ 
fer  employs  the  Korobkin  expression  for  the  Stanton  number,  Eq.  49.  The 
procedure  for  obtaining  the  heat  flux  from  this  expression  is  as  follows: 

2.  Compute  the  pressure,  Pe,  at  the  edge  of  the  boundary  layer  at 
the  point  0  from  Eq.  12  for  the  pressure  coefficient.  The  Newtonian 
pressure  distribution  may  be  used  to  advantage  over  the  surface  of  the 
hemisphere. 

3.  Compute  the  static  temperature,  Te,  at  the  outer  edge  of  the 
boundary  layer,  assuming  isentropic  flow  from  the  stagnation  point,  from 
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4.  Compute  the  velocity,  ue(fps)  ,  of  the  flow  at  the  outer  edge  of 
the  boundary  layer  from  the  flight  velocity,  V^fps),  and  the  enthalpies, 
he(Btu/lbm)  and  h^CBtu/lbm) ,  corresponding  to  the  temperatures  Te  and  T^. 

u  -  /v  2  +  50,100  /h  -  h  \ 
e  00  \  »  e/ 

5.  Compute  the  Reynolds  number  corresponding  to  the  local  flow 
conditions  at  the  outer  edge  of  the  boundary  layer  at  the  distance,  s(ft), 


Re 

s 


p  u  s 
e  e 


6.  If  Res  <  10^,  assume  laminar  boundary  layer.  If  Res  >  10^,  assume 
turbulent  boundary  layer. 

* 

7.  Compute  a,  the  recovery  factor,  r  ;  b,  the  adiabatic  wall  enthalpy; 
and  c,  the  reference  temperature,  T*(°R)  : 


a.  Recovery  factor 


* 


r 


f  *1/2 
Pr 


*1/3 

Pr 


for  laminar  flow 


for  turbulent  flow 


b.  Adiabatic  wall  enthalpy 


*  u 

*  e 

h  +  r - - — 

e  50,100 


Btu 

lbw 


c.  Reference  temperature 

T*  -  T  +  0.50  /T  -  T  \  +  0.22  /t  -  T  \ 

e  ^  w  e)  \  aw  ej 


where  the  adiabatic  recovery  temperature,  T^,  corresponds  to  the 
adiabatic  wall  enthalpy. 
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8.  Compute  the  density,  p*(lbm/fts),  and  viscosity,  y*(lbf-sec/ft2) , 
corresponding  to  the  local  pressure,  Pe(lbf/ft2)  and  reference  temperature, 
T*(‘R) . 

9.  Compute  the  Reynolds  number  based  on  reference  conditions  and 
the  running  length  along  the  surface  from  the  stagnation  point  to  the 
given  station: 


* 

Re 

s 


* 

P  Ufi8 

* 

V 


10.  For  laminar  flow  on  the  hemisphere 
equation,  Eq.  49,  for  the  Stanton  number  on 


surface,  use  Korobkin's 
a  sphere: 


Si 


/  *\-0*6 
0.763  (Prj 


11.  The  heat  flux  at  the  station  of  interest  Is  then  computed  from 

(h»  -  h») 

12.  For  turbulent  flow,  the  Eckert  flat -plate  equation  may  be  used 
to  compute  the  Stanton  number: 


Si 


0.0296 


for  105  <  Re  <  107  • 

—  s  — 


13.  When  a  wider  range  of  Reynolds  number,  Res,  is  encountered  fo'. 
turbulent  flow,  the  Schultz-Grunow  relation  may  be  used: 
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Heat  Transfer  in  the  Transitional  Region 

The  procedure  for  the  determination  of  the  heat  transfer  in  the 
transitional  region  between  the  laminar  and  turbulent  heat-transfer 
regions  is  outlined  as  follows: 

1.  The  Reynolds  number,  Reg,  based  upon  the  momentum  thickness 
may  be  chosen  as  the  parameter  to  represent  transition.  This  momentum 
Reynolds  number,  Reg,  based  on  local  fluid  properties  at  the  outer  edge 
of  the  boundary  layer  and  the  laminar  boundary  layer  momentum  thickness, 
may  be  computed  from  Ref.  57: 


0.22 

Cf/2 


0.22 


Transition  Is  assumed  to  occur  at  Reg  -  250.  Highly  polished  surfaces 
will  delay  transition;  hence,  higher  values  of  Reg  for  transition  will 
occur.  These  values  will  have  to  be  determined  experimentally. 


2.  The  constant  B  is  evaluated  from  the  laminar  flow  such  that  the 
local  friction  coefficient,  Cf/2,  remains  continuous  from  the  laminar  to 
turbulent  transition  region: 


at  transition 


3.  Calculation  of  the  skin-friction  and  the  momentum  Reynolds 
number  through  the  transition  region  proceeds  through  the  numerical 
solution  of  the  following  equations: 
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4.  The  heat-transfer  distribution  is  then  calculated  from 


«.  ‘  (Cf/2)  (Pre)'2/3  »,"e  (h„»  ‘  »„) 


5.  Finally,  the  local  NusselL  number  is  calculated  from 


Nu  -  Pr 


1/3  /  ^av 

h  -  h 
i  w 


(p./pQ(V\) 

(\) 


where 


Re 


p  u  R 
e 


is  the  Reynolds  number  based  on  flight  velocity,  stagnation  properties, 
and  the  radius  of  the  hemisphere. 


CONCLUSIONS 


PRESSURE  DISTRIBUTIONS 

In  addition  to  a  literature  survey,  methods  for  determining  the 
pressure  distribution  on  bodies  of  revolution  in  flight  regimes  from 
high  subsonic  through  supersonic  are  discussed  in  this  report  as  stated 
in  the  Introduction.  From  this  investigation,  the  modified  Newtonian 
and  the  modified  Newtonian/shock-expansion  methods  give  good  accuracy 
with  a  minimum  of  effort  for  the  following  cases: 

1.  For  a  hemisphere-cylinder,  the  pressure  distribution  may  bu 
determined  from  either  Eq.  3  or  13  in  combination  with  Eq.  12  for  the 
Mach  number  range  from  0.75  to  2.0,  and  from  Eq.  2  and  3  in  combination 
with  the  shock-expansion  method  for  Mach  numbers  above  2. 

2.  For  other  blunt  bodies  of  revolution  that  have  continuous 
curvature  in  the  nose  region,  Eq.  15  is  used  for  Mach  numbers  from 
0.7  to  2.0,  and  Eq.  2  is  used  for  Mach  numbers  above  2. 
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3.  For  polcted  bodies,  Eq.  6  and  the  modified  Newtonian/shock- 
expansion  method  are  used  for  Mach  numbers  above  1.0. 

4.  Equation  11  may  be  used  for  cylindrical  afterbodies. 

For  pointed  bodies  at  Mach  numbers  other  than  those  specified  as 
being  ammenable  to  the  modified  Newtonian/shock-expansion  method,  the 
most  useful  techniques  are  those  of  the  small-perturbation  theory  as 
developed  by  Van  Dyke  and  Spreiter. 


HEAT-TRANSFER  DISTRIBUTIONS 

This  report  gives  an  evaluation  of  applicable  blunt-body  heat- 
transfer  rates  for  high-subsonic  to  low-supersonic  flight.  An  outline 
of  an  engineering  method  of  analysis  for  the  prediction  of  blunt-body 
heat  transfer  within  this  flight  regime  is  also  presented.  From  this 
study,  the  following  conclusions  can  be  made. 

1.  Stagnation-point  heat  transfer  to  blunt-body  vehicles  for  a 
range  of  flight  Mach  numbers  from  2.3  to  3.1  has  been  predicted  by  the 
formula  of  Fay  and  Riddell  with  an  accuracy  from  8  to  20Z,  inclusively. 
However,  the  accuracy  of  prediction  of  wind  tunnel  data  by  the  formula 
of  Fay  and  Riddell  is  dependent  upon  the  enthalpy  difference  across  the 
boundary  layer  and  not  upon  the  tunnel  Mach  number.  Thus,  experimental 
datfi  are  required  to  determine  the  applicability  of  the  formula  cf  Fay 
and  Riddell  for  Mach  number  regions  of  less  than  2. 

2.  The  laminar  heat-transfer  distributions  over  blunt  bodies  of 
revolution  at  Mach  nusbers  between  approximately  2  and  5  are  best 
represented  by  Lees'  method.  If  laminar  flow  persists  to  the  juncture 
of  hemisphere  and  cone  on  a  hemispherically  capped  cone,  Lees'  method 
will  predict  the  effect  of  pressure  gradient  on  the  reduction  of  the 
heat  transfer.  However,  for  the  laminar  region  between  the  stagnation 
point  and  the  body  station  at  45  degrees,  little  difference  is  noted 
between  Korobkin's  method  and  Lees'  method. 

3.  The  method  of  predicting  the  heat-transfer  rate  in  the  transition 
region  on  blunt-body  vehicles  developed  by  Crasci,  Mackenzie,  and  Lirby; 
Nardo  and  Sadler;  and  Stetson  is  apparently  a  rational  and  fairly  accurate 
method  for  Including  transitional  behavior  in  heat-transfer  calcu  .atic?ns . 
Nevertheless,  methods  for  predicting  transition  still  reiuain  to  be 
developed. 

4.  Fully  developed,  turbulent  heat-transfer  rates  over  blunt-body 
vehicles  for  a  Mach  number  range  from  1.6  to  above  5  are  adequately 
represented  by  Eckert's  flat-plate  reference  enthalpy  technique.  No 
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turbulent  heat-transfer  data  are  available  In  the  open  literature  for 
either  flat  plates  or  for  blunt-body  systems  in  the  Mach  number  range 
from  0.7  to  1.6.  Such  data  will  be  required  to  provide  confidence  in 
the  application  of  the  engineering  methods  discussed  for  this  range  of 
Mach  numbers.  Extrapolation  of  the  flat-plate  formulas  into  this  range 
of  Mach  numbers  may  be  of  questionable  accuracy. 
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Appendix 

BIBLIOGRAPHY 

This  bibliography  consists  of  three  sections.  The  first  section 
is  an  annotated  bibliography  that  cites  and  critically  reviews  51  entries 
of  particular  interest  in  the  determination  of  pressure  and  heat-transfer 
distributions  on  bodies  of  revolution  for  subsonic,  transonic,  and 
supersonic  flight  regimes.  The  second  section,  a  general  bibliography, 
cites  516  entries  of  interest  in  this  area.  Some  additional  references 
covering  hypersonic  flow,  boundary-layer  theory,  real  gas  effects,  and 
corrections  for  angle  of  attack  or  yaw,  which  contain  pertinent  informa¬ 
tion,  are  also  included.  This  overall  bibliography  is  not  exhaustive, 
but  it  is  extensive  and  does  present  the  major  portion  of  those  articles 
having  a  direct  bearing  in  the  development  of  engineering  methods  of 
computing  pressure  and  heat-transfer  distribution. 

The  annotated  bibliography  contains  entries  arranged  alphabetically 
by  author.  Those  entries  covering  pressure  distribution  are  presented 
first  (Ref.  1—23),  followed  by  those  covering  heat-transfer  distribution 
(Ref.  24—51).  The  general  bibliography  contains  entries,  again,  arranged 
alphabetically  by  author;  however,  these  are  separated  into  six  main  groups 
by  subject.  The  largest  groups,  those  covering  pressure  distribution  and 
heat-transfer  distribution,  are  further  subdivided  for  more  convenient 
use. 
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ANNOTATED  BIBLIOGRAPHY 


PRESSURE  DISTRIBUTION 

1.  Adams,  M.  C.  and  Sears,  W.  R. ,  "Slender-Body  Theory — Extensions  and  Review," 
Journal  of  Aeronautical  Sciences,  Vol.  20,  No.  2,  February  1953,  pp.  85-98. 
(see  also,  Sears,  W.  R. ,  "Small  Perturbation  Theory"  High  Speed  Aerodynamic 
and  Jet  Propulsion,  Vol.  6,  Part  C. ,  Princeton  University  Press,  Princeton, 
New  Jersey,  1954. 

A  review  is  presented  of  small  perturbation  theory  as  applied  to  wings  and 
slender  bodies  of  revolution.  Cons'1  deration  is  given  to  determination  of 
pressure  distribution  from  theoretical  developments  which  are  obtained  by 
linearized  methods,  similarity,  end  simple  corrections  to  incompressible 
distributions.  Tills  report  provides  a  good  introduction  to  the  subject 
and  suggests  several  useful  modifications  of  earlier  methods  which  give 
improved  accuracy.  A  number  of  comparisons  of  theoretical  and  test  results 
are  given.  The  test  data,  however,  has  been  drawn  from  other  sources. 

2.  Andrews,  J.  S.,  "Steaay  State  Airload  Distribution  on  a  Hammerhead-Shaped 
Payload  of  a  Multistaged  Vehicle  at  Transonic  Speeds,"  D2-22947-1,  The 
Boeing  Company,  Seattle,  Washington,  February  1964. 

Data  is  presented  for  pressure  distribution  on  a  hemisphere  cylinder  for 
a  Mach  number  range  from  0.8  to  2.0  (see  Table  1A  for  data  summary).  Ad¬ 
ditional  and  empirical  modifications  of  the  modified  Newtonian  theory  is 
developed  which  correlates  the  data  for  0.75  <  M  <  2.0. 

C  -  C  sin2 6  +  F(0.,M  ) 
p  p  b  b  °» 

max 

F(0fa,Mj  “  [0-78(^)2’27  sin',0b  -  O.95e"2,235<VDlcos0b 
00 

3.  Baer,  A.  L. ,  "Pressure  Distribution  on  a  Hemisphere  Cylinder  at  Supersonic 
Mach  Numbers,"  AEDC-TN-61-96,  August  1961  I’AD  261  501). 

Experimental  pressure  distributions  were  obtained  over  the  hemispheric*! 
nose  of  AGARD  Model  E  at  Mach  numbers  ranging  from  2  to  8  (see  Table  1A 
for  summary  of  data).  Tests  results  are  compared  with  modified  Newtonian 
theory  and  show  good  agreement.  Comparison  with  the  numerical  methods  of 
Van  Dyke  and  Gordon  (NASA  TR  R-l,  1959)  are  also  given. 

4.  Blick,  E.  F.  and  Francis,  J.  E. ,  "Spherically  Blunted  Cone  Pressure  Distri¬ 
bution,"  AlAA  Journal,  Vol.  4,  No.  3,  March  1966,  pp.  547-549. 

The  modified  Newtonian-shock  expansion  method  is  applied  to  the  determina¬ 
tion  of  the  pressure  distribution  on  spherically  blunted  cones  The  angle 
at  which  the  modified  Newtonian  and  the  shock-expansion  methods  match  is 
given  by  the  empirical  equation 
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6  -  0.349  +  [0.198  +  0.80(M  -  2.8)_0,4551 

ID  op  * 

arc  coo  [1.05/lnM^) ]rad 
for  the  range  3  <_  <_  15. 

5.  Cole,  J.  D. ,  "Newtonian  Flow  Theory  for  Slender  Bodies,"  Journal  of  Aero¬ 
nautical  Sciences,  Vol.  24,  No.  8,  June  1957,  pp.  448-455. 

Newtonian  flow  theory  is  examined  from  the  point  of  view  of  gas  dynamics 
and  hypersonic  small-disturbance  theory.  A  general  solution  of  the  first 
approximation  for  the  flow  past  slender  bodies  at  zero  angle  of  attack 
is  given. 

This  paper  is  one  of  the  earlier  efforts  in  application  of  Newtonian 
theory  and  its  major  concern  is  the  hypersonic  flow  regime. 

6.  Guderley,  G.  and  Yoshlhara,  H.,  "Axial-Symmetric  Transonic  Flow  Patterns." 
USAF  Technical  Report  5797.  USAFAMC ,  Dayton,  Ohio,  1949. 

A  solution  for  the  flow  field  about  an  axlsymmetrlc  body  at  free  stream 
Mach  number  one  is  presented.  The  method  is  of  the  indirect  type,  work¬ 
ing  from  an  assumed  potential  flow  field  to  the  resulting  description  of 
a  nose  shape.  This  approach  appears  to  be  too  complicated  to  be  useful 
in  terms  of  an  "engineering"  solution  for  pressure  distribution  on  a  speci¬ 
fic  body. 

7.  Katz,  J.  R. ,  "Pressure  and  Wave  Drag  Coefficients  for  Hemisphere,  Hemisphere 
Cones,  and  Hemisphere  Ogives,"  NAVORD  Report  5849,  March  1958. 

Experimental  pressure  coefficients  for  hemisphere  cylinders  and  hem*  sphere 
cones  are  presented  end  compared  with  modified  Newtonian  theory  (a  summary 
of  the  data  is  presented  in  Table  1A).  No  data  is  presented  for  pressure 
distribution  on  hemisphere  ogives.  The  results  show  good  correlation  of 
experimental  data  for  a  Mach  number  range  2.5  to  6.5  with  the  modified  New¬ 
tonian  form 

C  ■  C  sin20. 

P  P  b 

r  r  mflv 


8.  Loitone,  E.  V. ,  "Linearized  Subsonic  and  Supersonic  Flow  about  Inclined 
Slender  Bodies  of  Revolution,"  Journal  of  Aeronautical  Sciences,  Vol.  14, 
No.  11,  November  1947,  pp.  631-642. 

Rigorous  first-order  linearized  solutions  are  obtained  for  subsonic  and 
supersonic  flow  about  bodies  of  revolution.  The  surface  pressure  is  found 
by  means  of  a  Taylor's  Series  Expansion  in  terms  of  the  cross-sectional 
area.  The  resulting  expressions  are  in  satisfactory  agreement  with  exist¬ 
ing  experimental  data.  The  first  order  solution  for  the  subsonic  regime 
is  strictly  valid  only  for  incompressible  flow  or  for  a  slender  pointed 
body  in  the  compressible  subsonic  regime.  However,  except  in  the  region 
of  a  stagnation  point,  the  solution  is  also  a  valid  approximation  for 
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slender  bodies  with  blunt  ends.  The  solution  for  the  supersonic  regime  is 
only  valid  where  the  entire  flow  field  is  supersonic.*  This  imposes  the  re¬ 
quirement  that  the  nose  be  pointed  and  the  shock  wave  attached. 

Some  comment  on  the  pressure  distribution  at  small  angles  of  atta'dc  is  also 
made. 

9.  Laitcne,  E.  V.,  "The  Subsonic  Flaw  about  a  Body  of  Revolution,"  Quarterly 
Applied  Mathematics,  Vol.  5,  Part  2,  July  1947. 

A  solution  is  presented  for  the  linearized  flow  potential  around  axisynmetric 
bodies  of  revolution  in  the  subsonic  regime: 

4  +  4.  +  1/r  $  -  0 

rxx  Trr  yr 

This  is  accomplished  by  approximating  the  source-sink  distribution  integral 
by  a  Taylor  Series  Expansion.  This  process  yields  the  following  first-order 
approximation  for  the  pressure  distribution 

„  ,1s'  ,l-2x/d,  Br  ,  „  x,,  xv  , 


„  .Is  ,l-2x/dv  t. ,  Br  ,  „  x,,  xv  , 

W  -  •  a  + 1» - - 1-  2  ru  -  ) 

+  1-j-  Ld  -  ^  l2(i  -  T  +  +  •  •  •  •  -  <t’>2 

where  the  coordinates  are  as  follows 


The  values  of  C  obtained  are  strictly  valid  only  for  incompressible  flow 
or  for  pointed  ^slender  bodies  in  subsonic  compressible  flow.  However,  ex¬ 
cept  in  the  region  of  a  stagnation  point,  the  expression  offers  a  valid 
approximation  for  the  pressure  distribution  on  blunt  slender  bodies  of 
revolution. 

10.  Lees,  L. ,  "Hypersonic  Flow,"  Proceedings  of  the  Fifth  International  Aero¬ 
nautical  Conference,  IAS-RAES.  June  20-23,  1955. 

A  considerable  contribution  to  the  usefulness  of  Newtonian  theory  is  made 
by  suggesting  a  modified  form  for  the  pressure  coefficient  on  blunt  bodies 

C  ■  C  sinz8. 


where  C 


-  the  value  at  the  stagnation  point. 


11.  Love,  E.  S.,  "Generalized  Newtonian  Theory,"  Journal  of  Aerospace  Sciences. 
Vol.  26,  No.  5,  May  1959,  p.  314. 
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A  mean*  of  applying  modified  Newtonian  theory  to  pointed  bodies  of  revolution 
is  given.  Using  Lees'  form 


C 

P 


sin26 

D 


Cn  is  defined  as 
nnax 


C 

p 

max 


C  /sin2e 
PN 


N 


where 

C  *  C  at  nose  of  cone  with  semi-vertex  angle  8,, 

pN  p  » 

9  ,  ■  local  slope  at  pointed  nose  of  the  body 
N 

This  application  is,  however,  limited  to  the  portion  of  the  supersonic  flight 
regime  where  the  body  will  have  an  attached  shock  at  the  nose. 


12.  Matthews,  C.  W.,  "Comparison  of  the  Experimental  Subsonic  Pressure  Distribu¬ 
tion  about  Several  Bodies  of  Revolution  with  Pressure  Distribution  Computed 
by  Means  of  the  Linearised  Theory,"  NACA  TN  2519,  February  1952. 


Experimental  pressure  distributions  on  two  prolate  spheroids  of  fineness  rat¬ 
ios  6  and  10,  an  ogive  body  and  a  prolate  spheroid  with  an  angular  bump  are 
conparcd  with  theoretical  results  obtained  from  applying  the  Prandtl-Glauert 
correction  to  the  incompressible  potential  flow  equations  and  from  linearized 
theory.  It  is  concluded  that  the  rather  more  complicated  Prandtl-Glauert 
method  gives  satisfactory  results  except  in  the  region  immediately  adjacent 
to  the  ends.  It  is  also  concluded  that  the  linearized  theory  does  not  yield 
satisfactory  results.  However,  it  should  be  noted  that  in  arriving  at  this 
conclusion,  Matthews  used  an  incomplete  expression  for  Cp.  Later  in  review¬ 
ing  this  work  Adams  and  Sears  point  out  that  using  the  complete  expression 
for  Cp 


linearized  theory  predicts  the  experimental  results  of  this  report  equally 
as  well  as  the  Prandtl-Glauert  convection  method  and  with  less  effort. 

13.  Okauchi,  K. ,  "Summary  of  Results  of  Pressure  Coefficient  Calculations  for 
a  Spherical- Tipped  Tangent-Ogive  Body,"  NAVWEPS  Report  8048.  October  1962. 

A  review  of  results  obtained  in  work  conducted  by  the  Naval  Ordnance  Test 
Station  to  develop  procedures  for  calculating  pressure  on  practical  body 
shapes  is  presented. 

The  body  used  for  the  test  had  a  spherically-tipped  (rc  ■  0.094  in.)  tan¬ 
gent  ogive  nose  (r^  -  0.592  in.).  Comparison  of  test  results  and  theoreti¬ 
cal  results  obtained  from  Van  Dyke's  second-order  method  is  made  at  - 
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1.53  and  2.17.  In  Che  theoretical  analysis  it  is  assumed  that  the  body 
is  pointed.  The  agreement  is  quite  good  (away  from  the  nose  tip).  For 
M«  ■  3.24  and  4.84  a  method  of  characteristic  solution  is  compared  with 
test  results.  Again  the  body  is  assumed  to  be  pointed  in  the  theoretical 
calculations  and  the  results  showed  good  agreement. 

This  report  also  contained  data  for  Moo  ■  0.75,  0.85,  0.95,  1.05,  and  1.15. 
There  Is  no  attempt  to  present  theoretical  comparisons  for  these  Mach  num¬ 
bers. 

14.  Perry,  J.  C.  and  Pasiuk,  C.  ,  "A  Comparison  of  Solutions  to  a  Blunt  Body 
Problem,"  AIAA  Journal,  Vol.  4,  No.  8,  August  1966,  p.  1438. 

Consideration  is  given  to  the  various  numerical  methods  of  predicting  the 
flow  in  the  nose  region  of  a  blunt  body  moving  at  supersonic  speeds.  Four¬ 
teen  different  solutions  were  obtained  and  compared  witn  experimental  re¬ 
sults  for  an  elipsoide  of  revolution  (r  -  /2x  -  4xz)at  M*  ■  3.  This  compcxi- 
son  showed  several  numerical  methods  closely  predicted  the  measured  test 
results.  There  was  scatter  between  the  pressures  predicted  by  various  num¬ 
erical  solutions,  however.  This  scatter  was  considerable  in  the  nose  region. 

The  numerical  methods  used  are  based  on  one  or  the  other  of  three  different 
computational  approaches  developed  by  Belotserkorski ,  Van  Dyke,  and  Swigart. 

A  method  developed  at  NASA  by  Lomax  and  Inouge  using  Van  Dyke's  approach  and 
a  method  using  Belotserkorski' a  approach  developed  at  NWL  by  DIDonato  pro¬ 
vided  quite  good  results. 

15.  Sprelter,  J.  R. ,  "Aerodynamics  of  Wings  and  Bodies  at  Transonic  Speeds," 
Journal  of  the  Aerospace  Sciences,  Vol.  26,  No.  8,  August  1959,  pp.  465-486. 

A  brief  summary  is  presented  on  basic  concepts  and  the  principle  results 
that  have  emerged  from  numerous  experimental  and  theoretical  investigations 
of  transonic  flow  past  thin  wings  and  slender  bodies.  Emphasis  is  placed 
throughout  on  the  correlation  and  evaluation  of  resulted  provided  by  di¬ 
verse  methods  for  the  approximate  solutions  of  the  nonlinear  equations  of 
the  small  perturbation  theory  in  transonic  flow.  Also  a  critical  examina¬ 
tion  of  experimental  results  is  given. 

This  paper  is  particularly  valuable  as  a  review  of  material  pertinent  to 
transonic  flight.  However,  the  bulk  of  the  material  presented  is  for 
two-dimensional  bodies  or  pointed  bodies  of  revolution 

and  also  appears  in  Sprelter  and  Alksne  (NASA  TR  R-2)  (see  Table  la  for 
data  summary) . 

16.  Spreiter,  J.  R.  and  Alksne,  A.  Y.,  "Slender  Body  Theory  Based  on  Approxi¬ 
mate  Solution  of  the  Transonic  Flow  Equation,"  NASA  TR  R-2.  1959. 

Approximate  solutions  of  the  nonlinear  equations  of  the  small  perturbation 
theory  for  transonic  flow  are  found  for  pointed  slender  bodies  of  revolu¬ 
tion  for  flows  with  free-stream  Mach  number  1  and  for  flows  that  are  eith¬ 
er  purely  subsonic  or  supersonic.  These  results  are  obtained  by  applica¬ 
tion  of  a  method  based  on  local  linearization.  The  theory  is  developed 
for  pointed  bodies  of  arbitrary  shape  and  specific  results  are  given  for 
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cone  cylinders  and  parabolic  arc  bodies  at  zero  angle  of  attack.  Results 
are  compared  with  data  (see  Table  1A  for  data  summary). 

The  method  of  this  report  is  rather  complicated  but  does  provide  solutions 
for  the  pressure  distribution  on  pointed  bodies  in  transonic  flight  not 
found  elsewhere.  Additionally,  the  report  is  of  considerable  value  as  gen¬ 
eral  background  material. 

17.  Van  Dyke,  M.  D. ,  "The  Paraboloid  of  Revolution  Subsonic  Flow,"  Journal  of 
Mathematics  and  Physics.  Vol.  37,  No.  1,  January  1958,  pp.  38-51. 

The  solution  for  the  flow  about  a  paraboloid  of  revolution  is  developed 
and  it  Is  demonstrated  how  this  solution  may  be  applied  to  remedy  the  de¬ 
fects  which  arise  when  linearized  small  perturbation  theory  and  its  high¬ 
er  approximations  are  applied  to  bodies  having  round  noses.  This  report 
provides  useful  background  when  applying  the  methods  developed  by  Van  Dyke 
in  NASA  TR  R-47. 

18.  Van  Dyke,  M.  D. ,  "Practical  Calculation  of  Second-Order  Supersonic  Flow 
Past  Ncnlifting  Bodies  of  Revolution,"  NACA  TN  2744,  July  1952. 

Calculation  of  second-order  supersonic  flow  past  pointed  bodies  of  revolu¬ 
tion  at  zero  angle  of  attack  is  described  in  detail  and  reduced  to  routine 
computation.  Tables  of  basic  functions  and  standard  computation  forms  are 
presented.  The  procedure  is  summarized  so  that  one  can  apply  it  without 
necessarily  understanding  the  details  of  the  theory.  A  sample  calculation 
is  also  given. 

19.  Van  Dyke,  M.  D. ,  "Second-Order  Slender  Body  Theory — Axlsymmetrlc  Flow," 

NASA  TR  R-47.  1959. 

Small  perturbation  theory  for  subsonic  and  supersonic  flow  past  bodies  of 
revolution  is  developed.  Methods  are  developed  for  handling  the  difficul¬ 
ties  which  arise  at  round  ends  in  subsonic  flow.  In  addition  to  second- 
order  theory,  first  order  linearized  theory  and  the  transonic  approximation 
to  the  non-linear  equations  are  developed  and  discussed  in  some  detail. 

The  methods  presented  are  somewhat  difficult  in  application  but  the  report 
presents  a  very  good  comprehensive  review  of  perturbation  theory.  All 
results  obtained  are  compared  with  other  existing  theories  or  with  ex¬ 
perimental  data  (see  Table  1A  for  data  summary). 

20.  Van  Dyke,  M.  D. ,  "The  Similarity  Rules  for  Second-Order  Subsonic  and 
Supersonic  Flow,"  NACA  Report  1374,  1958. 

The  similarity  rules  for  linearized  compreosible  flow  theory  (Gothert's 
rule  and  its  supersonic  counterpart)  are  extended  to  second  order.  It 
is  shown  that  any  second-order  subsonic  flow  can  be  related  to  "nearly 
incompressible"  flow  past  the  same  body,  which  can  be  calculated  by  the 
Janzen-Rayleigh  method.  It  is  also  pointed  out  that  these  methods  are 
chiefly  useful  as  a  theoretical  aid  rcther  than  in  correlation  of  experi¬ 
mental  data.  This  arises  as  tests  on  four  related  bodies  would  be  needed 
to  isolate  the  four  functions  involved  in  the  similarity  expression. 
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21.  Vau  Dyke,  M.  D. ,  "A  Study  of  Hypersonic  Small  Disturbance  Theory,"  NACA  TO 
3173,  May  1954. 

Newtonian  theory  Is  suggested  as  a  practical  means  of  determining  the  pres* 
sure  distribution  over  bodies  of  revolution  in  hypereonlc  flight.  For  the 
present  purpose  the  value  of  this  report  Is  chiefly  historical  In  that  It 
records  the  first  application  of  Newtonian  theory  to  problems  of  current 
interest.  In  Its  original  form  Newtonian  theory  yields  a  pressure  coeffici¬ 
ent  of  the  form 

C  -  2  sin2e. 

P  b 

22.  Vendemia,  R.  J.,  "An  Engineering  Method  for  Rapid  Calculation  of  Supersonic- 
Hypersonic  Pressure  Distributions  on  Lifting  and  Non-Lifting  Pointed  Bodies 
of  Revolution  and  Several  Special  Cases  of  Blunt-Nosed  Bodies  of  Revolution," 
TG  752,  The  Johns  Hopkins  University  Applied  Physics  Laboratory,  Silver 
Spring,  Maryland,  November  1965. 

A  method  is  presented  for  calculating  pressure  distribution  and  local  Mach 
numbers  for  pointed  bodies  of  revolution  and  for  several  cases  of  blunt- 
nosed  bodies  in  the  supersonic  and  hypersonic  flight  regimes.  The  method 
used  is  a  hybrid  solution  involving  modified  Newtonian  and  shock-expansion 
theory.  This  approach  is  very  useful  and  provides  accurate  results  for  a 
variety  of  nose  shapes  and  fineness  ratloe  in  the  supersonic  to  hypersonic 
regime.  A  considerable  amount  of  test  duta  is  also  presented  and  compared 
with  theory  (see  Table  1A  for  a  summary  of  Mach  nunbers  and  nose  shapes). 

In  the  context  of  the  present  problem  this  report  is  of  considerable  value. 

23.  Wagner,  R.  D. ,  Jr. ,  "Some  Aspects  of  Modified  Newtonian  and  Prandtl-Meyer 
Expansion  Method  for  Axisymmetric  Blunt  Bodies  at  Zero  Angle  of  Attack," 
Journal  of  Aerospace  Scienrts,  Vol.  26,  No.  12,  Decenber  1959,  p.  851. 

A  method  is  presented  for  determining  the  point  at  which  the  modified  New¬ 
tonian  and  shock-expansion  solutions  should  be  matched.  0£p  and  are 
given  in  graphical  form  as  a  function  of  Mach  number  (3  <_  M  40)  for  a 
hemisphere  cylinder  in  air  and  helium. 


HEAT-TRANSFER  DISTRIBUTION 

24.  Beckwith,  I.  V.  and  Gallagher,  J.  J.,  "Heat  Transfer  and  Recovery 

Temperatures  on  a  Sphere  with  Laminar,  Transitional,  and  Turbulent  Bound¬ 
ary  Layers  at  Mach  Numbers  of  2.00  and  4.15,"  NACA  TN  4125.  December  1957. 

An  investigation  was  made  of  the  pressure,  equilibrium,  temperatures,  and 
aerodynamic  heat  transfer  on  a  sphere  and  on  the  nose  of  a  hemisphere  cylin¬ 
der  at  Mach  numbers  of  2.00  and  4.15.  The  heat-transfer  coefficients  at 
the  stagnation  point  from  the  tests  at  both  Mach  numbers  were  compared  with 
the  theory  of  Reshotko  and  Cohen  (NACA  TN  3513)  and  were  in  good  agreement. 
The  velocity  gradients  used  in  the  theory  were  evaluated  from  the  experi¬ 
mental  pressure  date. 
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The  local  heat-transfer  coefficients  from  the  tests  at  Mach  number  2.00 
were  In  agreement  with  the  theory  for  laminar  flow  in  the  region  up  to 
approximately  20*  from  the  stagnation  point.  Large  increases  in  heat- 
transfer  coefficient  caused  by  transition  to  turbulent  flow  occurred  beyond 
this  region  with  the  maximum  values  observed  at  about  40*  from  the  stagna¬ 
tion  point.  From  this  location  back  to  the  separation  point  the  approximate 
level  in  the  heat-transfer  coefficient  as  well  as  the  variation  with  distance 
was  fairly  well  predicted  by  a  simple  theory  based  on  turbulent  heat-transfer 
formulas  for  flat  plates. 

The  heat  transfer  coefficients  from  the  tests  at  Mach  number  4.15  at  the 
highest  Reynolds  numbers  were  genu-ally  similar  to  the  results  at  Mach  num¬ 
ber  2.00  except  that  transition  taually  occurred  farther  back  on  the  nose. 

The  maximum  values  of  the  heat  transfer  coefficient  downstream  of  transition 
were  again  in  reasonably  good  agreement  with  the  simple  turbulent-flow  theory. 
For  Reynolds  number  of  about  3.00  x  10&  or  less,  the  data  over  the  entire 
hemispherical  nose  were  in  agreement  with  the  theory  for  laminar  flow. 

25.  Brindle,  C.  C.  and  Malia,  M.  J.,  Jr.,  ’’Longitudinal  Aerodynamic  and  Heat 
Transfer  Characteristics  of  a  Hemisphere-Cylinder  Missile  Configuration 
with  an  Aerodynamic  Spike,"  Aero  Report  No.  1061.  Washington,  D.C.,  The 
David  Taylor  Model  Basin  Aerodynsmlcs  Laboratory,  AD  No.  417  466,  July 
1963. 

A  hemisphere  cylinder  with  an  aerodynamic  spika  (produced  by  a  small  nos- 
xle  exhausting  air  from  the  hemispherical  nose)  was  tested  in  the  18-lnch 
wind  tunnel  at  Mach  nuabars  of  1.87  and  2.16.  Heat-transfer  coefficients 
were  measured  for  both  hot  and  cold  bodies.  No  comparison  or  correlations 
with  theory  were  obtained 

26.  Brunk,  J.  E.  and  White,  W.  N. ,  "Free-Flight  Hypersonic  Heat  Transfer  and 
Boundary  Layer  Transition  Studies:  HTV  Flights  A-40  and  A-41,”  ARL  Tech¬ 
nical  Report  60-311,  Wright-Patterson  Air  Force  Base,  Ohio,  Aeronautical 
Research  Laboratory,  December  1960. 

Two  HTV-1  hypersonic  test  vehicles,  Rounds  A-40  and  A-41,  were  flown  at 
Holloman  AFB  in  October  1959  with  blunted  and  sharp  20*  half  angle  nose 
cones,  respectively.  A  maximum  flight  velocity  of  5800  feet  per  second 
was  attained,  corresponding  to  a  local  sharp  core  Mach  number  and  unit 
Reynolds  number  of  3.4  and  50  x  10®  per  foot  respectively. 

Surface  temperatures  and  aerodynamic  heating  rates  were  obtained  for  one 
station  and  three  radial  positions  on  the  conical  portion  of  the  blunted 
nose  cone  (Round  A-40)  end  at  three  stations  on  each  of  two  longitudinal 
rays  on  the  sharp  cone  (Round  A-41).  Uncertainties  in  the  data  reduction 
Invalidated  the  laminar  results  but  the  heating  rates  obtained  in  turbu¬ 
lent  regions  on  the  conical  sections  showed  good  agreement  with  the  theo¬ 
retical  heating  rates  based  upon  Eckert's  flat -plate  reference  enthalpy 
method.  The  surface  finish  was  specified  as  32  micrc  inches  (RMS)  maximum. 

27.  Buglia,  J.  J.,  "Heat  Transfer  and  Boundary-Layer  Transition  on  a  Highly- 
Polished  Hemisphere-Cone  in  Free  Flight  at  Mach  Numbers  up  to  3.14  and 
Reynolds  Numbers  up  to  24  x  106,"  NASA  TN  D-955,  September  1961. 
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The  model  used  for  the  flight-testing  had  a  ratio  of  nose  radius  to  base 
radius  of  0.74  and  a  half-angle  of  14.5*  and  was  flight-tested  at  Mach 
numbers  up  to  4.70.  The  nose  of  the  model  had  a  surface  roughness  of  2  to  5 
microinches  as  measured  with  an  interferometer.  A  comparison  with  results 
from  rougher  models  indicated  that  the  high  degree  of  polish  was  instrumental 
in  delaying  the  transition  from  laminar  to  turbulent  flow.  However,  data  were 
obtained  at  Mach  numbers  of  2.32,  2.47,  2.63,  2.80,  2.97,  and  3.14  and  prob¬ 
ably  represent  the  most  complete  set  yet  obtained  for  correlation  with  theo¬ 
retical  methods.  Stagnation  point  heating  values  were  in  best  agreement  with 
the  prediction  of  Cohen  while  Lees  theory  best  represented  the  distribution 
of  heat  transfer.  The  turbulent  flow  over  the  conical  portion  of  the  vehicle 
correlates  best  with  the  Eckert  flat-plate  reference  enthalpy  method. 

28.  Carter,  H.  S.  and  Bressette,  W.  E.,  "Heat  Transfer  and  Pressure  Distribution 
on  Six  Blunt  Noses  at  Mach  Number  of  2,"  NACA  RM  L57C18.  April  1957. 

The  heat  transfer  and  pressures  on  the  surfaces  of  six  blunt-nose  models  are 
presented  for  0°  and  5*  angles  of  attack.  The  tests  were  made  in  a  free  jet 
at  a  Mach  number  of  2  with  a  total  temperature  of  935*R  and  a  Reynolds  num¬ 
ber  per  foot  of  about  14  x  106.  The  authors  state  that  the  heat  transfer  over 
the  conical  portion  of  the  models  was  correlated  with  Van  Drlest's  flat-plate 
theory  while  the  heat  transfer  at  the  stagnation  point  was  compared  with  the 
theory  of  Reshotko  and  Cohen.  However,  the  data  does  not  compare  well  with 
the  results  of  Van  Driest  when  comparison  is  based  upon  local  conditions  at 
the  outer  edge  of  the  boundary  layer.  Insufficient  data  precluded  comparison 
on  the  reference-enthalpy  method. 

29.  Carter,  H.  S.  and  Carr,  R.  E. ,  "Free-Flight  Investigation  of  Heat  Transfer  to 
an  Unswept  Cylinder  Subjected  to  an  Incident  Shock  and  Flow  Interference  from 
an  Upstream  Body  at  Mach  Numbers  up  to  5.50,"  NASA  IN  D-988.  October  1961. 

Skin-tewperature  measurements  were  made  along  the  stagnation  line  of  the 
transverse  cylinder  of  a  cruciform  test  configuration  constructed  of  two 
nearly  equal  size  tubes.  The  flow  on  this  stagnation  line  was  influenced 
by  boundary- layer  separation  and  by  end  effects  as  well  as  by  the  intersec¬ 
ting  shocks.  Under  these  conditions  the  measured  heat-transfer  rates  had 
inconsistent  variations  both  in  magnitude  and  distribution  which  precluded 
correlation  of  the  heating  rates  with  position  of  the  intersecting  shocks 
or  with  theory. 

30.  Chauvin,  L.  T.,  "Pressure  Distribution  and  Pressure  Drag  for  a  Hemispheri¬ 
cal  Nose  at  Mach  Numbers  2.05,  2.54,  and  3.04,"  NACA  RM  L52K06.  December  1952. 

The  results  of  an  experimental  investigation  of  the  pressure  distributions  on 
a  hemispherical  nose  3.98  inches  in  diameter,  mounted  or.  a  cylindrical  support, 
were  presented  for  wind  tunnel  Mach  numbers  of  2.05,  2.54.  and  3.04  and  for 
Reynolds  numbers  of  4.44  x  106,  4.57  x  106,  and  4.16  x  10&,  respectively,  bas¬ 
ed  on  body  diameter  and  free-stream  conditions. 

This  report  presents  basic  tunnel  operating  parameters  required  for  the 
evaluation  of  the  convective  heat  transfer  results  presented  by  Chauvin 
and  Maloney  in  NACA  RM  L53L08a  (February  1954). 
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31.  Chauvin,  L.  T.  and  Maloney,  J.  P.t  "Experimental  Convective  Heat  Transfer 
to  a  4-Inch  and  6- Inch  Hemisphere  at  Mach  Number*  from  1.62  to  3.04," 

HACA  «h  L53L08a.  February  1954. 

Equilibrium  skin  temperatures  and  convective  heat  transfer  were  measured  for 
a  4-inch  and  4- Inch  diameter  hemispherical  nose  in  the  Langley  ei-inch  and 
12-inch  blowdown  wind  tunnels,  respectively.  The  Mach  number  f  ir  the  4-lnch 
nose  ranged  from  1.62  to  3.04  for  a  Reynolds  number  based  on  body  diameter 
of  approximately  4.5  x  10®.  The  6-inch  nose  vas  tested  at  M  -  1.99  at  a 
Reynolds  nunber  based  on  body  diameter  of  6.4  x  106  and  for  two  surface  con¬ 
ditions  . 

The  measured  Stanton  numbers  correlated  well  with  the  Korobkin  equation 
over  the  spherical  section  of  the  hemisphere  for  the  laminar  region.  In 
the  turbulent  region  the  flat-plate  reference  method  of  Eckert  represented 
the  data  fairly  well.  Transition  vas  observed  at  Reynolds  number  of  ap¬ 
proximately  1  x  106  based  on  distance  from  the  stagnation  point  and  condi¬ 
tions  Just  outside  the  boundary  layer  and  corresponding  to  a  regi  on  between 
the  45*  and  the  60*  stations. 

32.  Cohen,  N.  B.,  "Boundary-Layer  Similar  Solutions  and  Correlation  Equation 
for  Laminar  Keat-Transfer  Distribution  in  Equlpibrium  /dr  at  Velocities 
up  to  41,000  Feet  Per  Second,"  NASA  TI  R-118.  1961. 

Exact  solutions  for  flat-plate  flow,  axis  yarn*  trie  stagnation  flow  and  stag¬ 
nation  flow  for  a  yawed  infinite  cylinder  are  obtained  and  correlating  formu¬ 
las  for  the  appropriate  heat  transfer  functions  are  developed.  A  correla¬ 
tion  formula  for  the  heat  transfer  distribution  function  for  an  arbitrary 
favorable  pressure  gradient  on  a  body  of  revolution  or  yawed  infinite  cylin¬ 
der  is  determined  for  application  to  various  local  similarity  methods.  Per¬ 
tinent  equations  are  collected  in  a  single  section  for  convenience. 

The  methods  presented  in  this  report  should  be  extended  to  lower  Mach  num¬ 
bers.  The  correlation  presented  agrees  better  with  experimental  data  than 
other  well-known  formulas  and  the  range  of  application  should  be  extended. 

33.  Crawford,  0.  H.  and  McCauley,  W.  D. ,  "Investigation  of  the  Laminar  Aerody¬ 
namic  Heat-Transfer  Characteristics  of  a  Hemisphere  Cylinder  in  the  Langley 
11-Inch  Hypersonic  Tunnel  at  a  Mach  Number  of  6.8,"  NACA  Report  1323.  1957. 

Aerodynamic  heat  transfer  data  to  a  nonisothermol  hemisphere  cylinder  evalu¬ 
ated  in  the  Langley  11- inch  hypersonic  tunnel  at  a  Mach  number  of  6.8  and 
a  Reynolds  number  from  approximately  0.14  x  10&  to  1.06  x  lh^based  on  dia¬ 
meter  and  free-atream  conditions  were  reported.  The  experimental  heat-trans¬ 
fer  coefficients  were  slightly  less  over  the  whole  body  than  those  predicted 
by  the  theory  of  Stine  and  Hanlass  (NACA  TN  3344)  for  an  Isothermal  surface. 
For  stations  within  45*  of  the  stagnation  point  the  heat-transfer  coeffici¬ 
ents  could  be  correlated  by  a  single  relation  between  local  itanton  number 
and  local  Reynolds  number. 

34.  Creed,  R.  J.,  Mackenzie,  D.  A.  and  Libby,  P.  A.,  "An  Investigation  of 
laminar.  Transitional,  and  Turbulent  Heat  Transfer  on  Blunt-Nosed  Bodies 
in  Hypersonic  Flow,"  Journal  of  the  Aerospace  Sciences,  Vol.  27,  No.  6, 

June  1960,  pp.  401-414. 
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Experimental  results  for  1  seiner,  transitional,  end  turbulent  besting 
rates  measured  by  of  the  shrouded  model  technique  were  reported. 

The  Reynolds  number  was  varied  over  a  range  from  C.4  x  10®  to  4.0  x  10®  while 
the  stagnation  to  wall  enthalpy  ratio  from  2.3  to  approximately  1.5. 

The  experimental  data  were  compared  to  the  laminar  hypersonic  boundary  layer 
theory  of  Lees  and  shown  to  be  in  good  agreement  on  the  conical  portion  of 
the  model.  On  the  spherical  portion  of  the  model  the  data  are  approximately 
20  per  seat  higher  than  the  theoretical  prediction,  due  apparently  to  radia¬ 
tion  to  the  model. 

The  experimental  date  presented  in  this  study  is  not  directly  applicable  to 
the  problem  of  determining  the  aerodynamic  heat  transfer  to  blunt  bodies  at 
lover  supersonic  Mach  numbers.  However,  the  method  of  analysis  should  be 
explored  for  application  in  the  transitional  flow. 

Extensive  data  -re  also  presented  for  comparison  of  the  flat-plate  reference- 
enthalpy  theory  for  the  prediction  of  turbulent  heat  transfer  to  blunt  bodies. 
It  should  be  noted,  however,  that  these  date  were  obtained  for  stagnation 
temperatures  in  the  range  of  1500 *R  to  1700*R  which  correspond  tc  flight 
Mach  nunfcers  of  10  and  shove. 

35.  Dlaconis,  N.  S.,  Wisniewski,  R.  J.  ,  and  Jack,  J.  R. ,  "Heat  Transfer  and 
Boundary-Layer  Transition  on  Two  Blunt  Bodies  at  Mach  Humber  3.12,"  NACA  TN 
4099,  October  1957. 

Local  heat  transfer  values  were  measured  on  a  hemisphere  cone  cylinder  of 
diameter  of  1.4  laches  and  a  120*  Included  angle  cone  cylinder  at  a  Mach 
number  of  3.12.  Experimental  valuta  of  the  heat  transfer  distribution  for 
the  blunt  body  ..'greed  well  w.  tb  tne  heat  transfer  distribution  prediction 
of  the  theory  of  Lees  when  the  surface  roughness  was  less  than  16  micro- 
inches.  An  insufficient  amount  of  wall  temperature  information  prohibited 
comparison  of  the  distribution  results. 

36.  Hermann,  R.  and  Melnick,  W.  L. ,  "Aerodynamic  end  Heat  Transfer  Studies  with 
Evaporative  Film  Cooling  at  Hypersonic  Mach  Nunbers,"  RAL  KR  189,  University 
of  Minnesota,  Rosemount,  Minnesota,  Rosemount  Aeronautical  Laboratories, 
September  1962. 

In  order  to  establish  the  insulating  characteristics  of  an  evaporating  film, 
measurements  of  the  aerodynamic  heat  transfer  to  a  hemiepherically- capped 
cylinder  were  made  at  free-stream  Mach  numbers  of  0.30  and  6.8  without  mass 
injection.  The  low  speed  results  (M®  “  0.30)  indicated  that  the  heat  trans¬ 
fer  conductance  (h)  were  constant  In  the  stagnation  region,  in  agreement  with 
well-known  theories  (where  Nug  ■  0.93  The  laminar  heat  transfer 

distribution  on  the  hemisphere  cylinder  at  M®  ■  6.8  with  a  Revnr  Ids  number 
of  283,000  per  foot  was  well  represented  by  the  theory  nf  Kemp  Rose,  and 
Detra. 


37.  Kemp,  N.  H. ,  Rose,  P.  H.,  and  Detra,  R.  W.,  "Laminar  Heat  Transfer  Around 
Blunt  Bodies  in  Dissociated  Air,"  Journal  of  the  Aerospace  Sciences,  Vol. 
26,  No.  7,  July  1959,  pp.  421-430.' 
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An  extension  of  the  method  of  Fey  end  Riddell  le  presented  for  predicting 
laminar  heat- transfer  retes  to  blunt,  highly-cooled  bodies  with  constant 
well  tempereture  In  dlssocieted  elr  flow. 

Experimental  heat-transfer  rates  obtained  from  shock  tube  facilities  were 
presented  and  compared  with  the  theoretical  results.  Approximately  75  data 
points  were  presented  for  the  heat  transfer  distribution  over  a  hemispheri- 
cally-capped  cylinder.  These  data  were  obtained  In  a  shock  tube  at  ehock 
Mach  numbers  from  approximately  8  to  approximately  14,  and  hence  the  result 
lng  experimental  values  of  stagnation  point  heat  transfer  and  heat  transfer 
distribution  are  not  applicable  for  flight  Mach  numbers  less  than  4. 

38.  Korobkln,  I.,  "Local  Flow  Conditions,  Recovery  Factors  and  Heat-Transfer 
Coeffl  :ients  on  the  Nose  of  a  Hemisphere-Cylinder  at  Mach  Number  of  2.8," 
NAVORL  Report  286S,  White  Oak,  Maryland,  U.S.  Naval  Ordnance  Laboratory, 

May  1953. 

a 

Recovery  temperatures  and  heat  transfer  coefficients  were  measured  on  an 
one-inch  diameter  hemlsphe.rically  capped  cylinder  in  the  NOL  40  x  40  cm 
Aeroballistlco  Tunnel  No.  3  at  a  Mach  number  of  2.8.  The  locd.  heat-transfer 
coefficient  for  supersonic  flow  over  the  hemisphere  with  a  laminar  boundary 
layer  were  correlated  with  the  Incompressible  expression  presented  by  Sibulkln: 

N  -  K(8)  R  °'5  Pr0,A 
“D  *D 

39.  Korobkln,  I.,  "Laminar  Heat  Transfer  Characteristics  of  a  Hemisphere  for  a 
Mach  Number  Range  1.9  to  4.9,"  NAVORD  Report  3841,  White  Oak,  Maryland, 

U.S.  Naval  Ordnance  Laboratory,  October  1954. 

Pressure  and  heat  transfer  measurements  were  made  on  the  hemispherical  nose 
of  a  two-inch  diameter  hemisphere  cylinder  with  laminar  boundary  layer  flow. 

The  tests  were  conducted  in  the  NOL  40  x  40  cm  Aeroballistlcs  Tunnel  No.  i. 

Heat  transfer  measurements  were  made  In  the  Mach  number  ‘:£'ige  from  1.90  to 
4.87.  The  resultant  Stanton  number  correlates  with  the  expression 
for  the  stagnation  region  of  a  sphere 

CH  -  0.763  (Pr)"°-6  (Res)'0,5 

where  Pr  *  0.70  and  the  Reynolds  number  is  based  upon  the  local  conditions 
at  the  outer  edge  of  the  boundary  layer  and  the  length  parameter  is  the 
wetted  surface  distance  from  the  stagnation  point.  These  data  are  presen¬ 
ted  in  Figure  19. 

40.  Marvin,  J.  G.  and  Deiwert,  G.  S.,  "Convective  Heat  Transfer  in  Planetary 
Gases,"  NASA  TR  R-224,  July  1965. 

F.quilibrium  convective  heat  transfer  in  several  real  gases  was  investigated. 

The  gLS  considered  were  air,  nitrogen,  hydrogen,  carbon  dioxide,  and  argon. 
Solutions  ti  the  similar  form  of  the  boundary- layer  equations  were  obtained 
for  flight  velocities  from  10,000  to  30,000  ft/sec  for  a  range  of  parameters 
sufficient  to  define  the  effects  of  pressure  level,  pressure  gradient, 
boundary- layer  edge  velocity  and  wall  temperature.  Results  are  presented 
for  stagnation  point  heating  and  for  heating-rate  distribution. 
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Although  the  range  of  velocities  considered  in  this  reference  are  beyond 
those  of  this  bibliography,  the  method  of  analysis  and  trends  Indicated 
may  be  useful  in  developing  methods  of  numerical  solution  of  lower  Mach  num¬ 
ber  heat  transfer  problems. 

41.  Nardo,  S.  V.  and  Sadler,  R.  W. ,  "Heat  Transfer  and  Temperature  Distribution 
in  a  Hemispherical  Nose  Cone  in  Hypersonic  Flow,"  PIBAL  Report  No.  597,  Poly¬ 
technic  Institute  of  Brooklyn,  New  York,  Brooklyn,  New  York,  January  1962. 

The  heat  transfer  to  a  hemispherical  nose  cone  subjected  to  hypersonic  flow 
conditions  was  calculated  by  the  theory  of  Lees  for  Reynolds  numbers  from 
one  to  five  million.  Comparison  with  results  of  test  runs  in  the  Polytech¬ 
nic  hypersonic  tunnel  facility  Indicated  good  agreement.  The  experimental 
procedure  employed  the  shrouded  model  technique  which  is  described  in  the 
paper  of  Cresci,  Mackenzie,  and  Libby  (Reference  34). 

This  report  presents  further  details  of  a  computational  procedure  for  deter¬ 
mining  the  heat  transfer  in  the  transition  region  from  laminar  to  turbulent 
boundary  layers  on  a  hemispherical  nose.  Further  experimental  verification 
of  the  procedure  is  also  presented. 

42.  Neal,  L. ,  Jr.  and  Bertram,  M.  H.,  "Turbulent  Skin  Friction  and  Heat  Transfer 
Charts  Adapted  from  the  Spalding  and  Chi  Method,"  NASA  TN  D-3969,  May  1967. 

Charts  are  presented  which  allow  a  rapid  determination  of  local  and  average 
skin  friction  and  heat  transfer  on  flat  plates  in  air.  The  charts  cover  a 
Mach  number  range  from  0  to  20,  a  Reynolds  number  range  from  1()5  to  10^  in 
decade  increments  of  the  exponent  and  a  wall  temperature  stagnation  temp¬ 
erature  ratio  range  of  0.1  to  the  adiabatic  wall  case. 

The  accuracy  of  the  charts  has  not  been  established  in  the  Mach  number  range 
from  0.8  to  5  since  the  experimental  heat  transfer  data  used  to  support  the 
method  employed  in  the  development  of  the  charts  was  obtained  in  the  Mach 
number  range  from  about  5  to  9. 

43.  Sands,  N.  ar.d  Jack,  J.  R. ,  "Preliminary  Heat  Transfer  Studies  on  Two 
Bodies  of  Revolution  at  Angle  of  Attack  at  a  Mach  Number  of  3.12,"  NACA  TN 
4378,  September  1958. 

The  axial  temperature  and  Stanton  number  distributions  for  two  pointed  bodies 
of  revolution  at  zero  angle  of  attack  are  presented.  No  comparisons  were 
made  to  theory. 

44.  Schumann,  E.  E.  H. ,  "Engineering  Methods  for  the  Analysis  of  Aerodynamic 
Heating,"  RAD-TM-63-68,  Wilmington,  Massachusetts,  Research  and  Advanced 
Development  Division,  Avco  Corporation,  November  1963. 

Various  methods  cf  analysis  of  the  convective  heating  to  a  body  are  review¬ 
ed  for  flight  Mach  numbers  above  5.  This  report  is  of  value  in  comparing 
methods  of  analysis  with  experiment  for  Mach  numbers  from  8.5  to  16  but 
does  not  contribute  to  evaluation  of  heating  effects  below  Mach  numbers  of 
5. 
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45.  Sei£f,  A. ,  "Examination  of  the  Existing  Data  on  the  Heat  Transfer  of  Turbu¬ 
lent  Boundary  Layers  at  Supersonic  Speeds  from  the  Point  of  View  of  Reynolds 
Analogy,"  NACA  TN  3784,  August  1954. 

Heat  transfer  data  from  four  wind  tunnel  experiments  and  two  free-flight  ex¬ 
periments  with  turbulent  boundary  layers  were  compared  with  Reynolds  analogy 
and  agreement  was  found  within  eight  per  cent.  The  data  covered  the  Mach 
number  range  from  1.4  to  3.2  at  effective  Reynolds  numbers  from  0.4  million 
to  24  million.  They  were  obtained  on  a  variety  of  body  shapes  including 
flat  plates,  cones,  and  pointed  slender  bodies  of  revolution.  Extremely 
blunt  shapes  such  as  spheres  were  excluded  from  the  correlation. 

46.  Sommer,  S.  C.  and  Short,  B.  J.,  "Free  Flight  Measurements  of  Turbulent- 
Boundary  Layer  Skin  Friction  in  the  Presence  of  Severe  Aerodynamic  Heating 
at  Mach  Numbers  from  2.8  to  7.0,"  NACA  TN  3391.  1955. 

Experimental  measurements  of  average  skin  friction  of  the  turbulent  bound¬ 
ary  layer  were  presented  for  a  free-flying,  hollow-cylinder  model  at  Mach 
numbers  of  2.8,  3.8,  5.6,  and  7.0  at  conditions  of  high  rates  of  heat 
transfer.  No  heat  transfer  data  or  correlations  are  presented. 

47.  Stine,  H.  A.  and  Uanlass,  K. ,  "Theoretical  and  Experimental  Investigation 
of  Aerodynamic  Heating  and  Isothermal  Heat  Transfer  Parameters  on  a  Hemi¬ 
spherical  Nose  with  Laminar  Boundary  Layer  at  Supersonic  Mach  Numbers," 

NACA  TN  3344.  December  1954. 

Measurements  of  the  heat-transfer  parameter  Nu/*^  baaed  on  local  flow  con¬ 
ditions  just  outside  the  boundary  layer  and  length  of  boundary- layer  run 
were  carried  out  on  the  heated,  isothermal  surface  of  a  hemisphere  cylin¬ 
der  for  a  Mach  number  of  1.97  and  Reynolds  number  based  on  body  diameter 
from  0.60  x  10^  to  2.28  x  10^.  An  approximate  method  for  calculating  the 
distribution  of  the  heat-transfer  parameter,  applicable  to  any  body  of 
revolution  was  developed  which  predicts  within  about  +  18  per  cent  the 
experimental  results.  Since  the  data  presented  in  the  above  study  were 
obtained  from  an  internally  heated  body,  they  were  not  subjected  to  an¬ 
alysis  and  correlation  in  our  study. 

48.  Swanson,  A.  G.,  Buglla,  J.  J.,  and  Chauvin,  L.  T. ,  "Flight  Measurements 
of  Boundary-Layer  Temperature  Profiles  on  a  Body  of  Revolution  (NACA  RM- 
10)  at  Mach  Numbers  from  1.2  to  3.5,"  NACA  TN  4061,  July  1957. 

A  parabolic  body  of  revolution  with  a  pointed  nose,  the  NAG K  RM-10  was 
flight  tested  at  Mach  numbers  up  to  3.5  and  Reynolds  numbeis  up  to  140  x 
106  (based  on  length  to  the  probe  and  rake  measuring  station).  Skin- 
friction  coefficients  determined  by  the  momentum  method  from  the  tctai- 
pressure  rake  data  were  in  fair  agreement  with  the  Van  Driest  theory  for 
flat  plates  with  compressible  turbulent  boundary  layers.  Heat  transfer 
measurements  indicated  a  transition  Reynolds  number  of  about  15  x  10&. 

49.  Wilson,  R.  E.,  "Viscosity  and  Heat  Transfer  Effects,"  (Sections  13  and 
14,  Handbook  of  Supersonic  Aerodynamics),  NAVORD  Report  1488,  Vol.  5, 

August  1966. 
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This  report  presented  a  review  and  discussion  of  nearly  two  hundred  papers 
pertaining  to  the  study  of  viscous  effects  in  high-speed  flows.  The  empha¬ 
sis  in  this  volume  is  placed  on  the  presentation  of  theoretical  and  experi¬ 
mental  results  which  supposedly  can  be  used  by  the  missile  designer  for  the 
calculation  of  skin  friction  and  heat  transfer. 

The  method  of  calculation  suggested  in  this  report  for  the  heat  transfer 
to  blunt  bodies  of  revolution  is,  however,  incomplete  in  this  discussion 
of  stagnation  point  heat  transfer  and  the  distribution  of  heat  transfer 
over  the  body  surface.  Results  are  presented  for  high  Mach  numbers  but 
the  report  does  not  extend  the  analysis  to  Mach  numbers  less  than  5  but 
greater  than  0.87. 

50.  Hinkler,  E.  M.  and  Daoberg,  J.  E. ,  "Heat  Transfer  Characteristics  of  a 
Hemisphere  Cylinder  at  Hypersonic  Mach  Numbers,"  NAVORD  Report  4259,  White 
Oak,  Maryland,  U.S.  Navul  Ordnance  Laboratory,  April  11,  1957. 

The  heat-transfer  characteristics  of  the  laminar  compressible  boundary  layer 
on  a  hemisphere  cylinder  were  presented  for  free-stream  Mach  numbers  of  5, 
6.5,  and  8.  These  data  have  not  been  included  in  this  bibliography.  How¬ 
ever,  it  should  be  noted  that  excellent  correlation  was  attained  with  Van 
Drlest's  stagnation  region  expression  for  incompressible  flow. 

51.  Zoby,  E.  V.  and  Sullivan,  E.  M. ,  "Correlation  of  Free-Fllght  Turbulent 
Heat  Transfer  Data  from  Axlsymmetric  Bodies  with  Compressible  Flat-Plate 
Relationships,"  NASA  TN  D-3802,  January  1967. 

Published  experimental  turbulent  heat  transfer  data  obtained  over  a  range 
of  free-flight  conditions  and  body  shapes  (free-stream  Mach  numbers  from 
2.9  to  13.4,  free-stream  Reynolds  numbers  per  foot  from  0.64  x  10^  to  30.7  x 
106,  body  shapes  from  a  hemisphere  cylinder  to  a  sphere  cone  with  a  half 
angle  of  25*)  were  compared  with  calculated  turbulent  flat-plate  values. 

The  calculated  values  were  evaluated  by  use  of  a  modified  Reynolds  analogy 
and  the  skin-friction  relationships  of  Blasius  or  Schultz-Grunow  with 
compressibility  effects  accounted  for  by  evaluating  the  flow  properties 
on  reference  conditions.  For  reference  Reynolds  numbers  less  than  107,  the 
calculated  heating  rates  based  on  either  of  the  two  methods  correlated  well 
with  the  experimental  data.  For  reference  Reynolds  numbers  greater  than  107 
and  less  than  8  x  10?,  the  calculated  heating  rates  billed  on  the  Schultz- 
Grunow  relation  compere  better  with  the  available  experimental  data. 

Zoby  and  Sullivan  point  out  that  available  flight-test  data  correlate  with¬ 
in  20  per  cent  of  the  correlation  line  whereas  wind  tunnel  data  spread  out 
to  +  50  per  cent  of  the  correlation  line.  This  point  should  be  pursued  in 
further  study. 
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